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The quantum dynamics of mesoscopic or macroscopic systems is always complicated by their 
coupling to many "environmental" modes. At low T these environmental effects are dominated 
by localised modes, such as nuclear and paramagnetic spins, and defects (which also dominate the 
entropy and specific heat). This environment, at low energies, maps onto a "spin bath" model. 
This contrasts with "oscillator bath" models (originated by Feynman and Vernon) which describe 
delocalised environmental modes such as electrons, phonons, photons, magnons, etc. The couplings 
to N spin bath modes are independent of N (rather than the ~ 0{l/\fN) dependence typical of 
oscillator baths), and often strong. One cannot in generalmap a spin bath to an oscillator bath (or 
vice- versa); they constitute distinct "universality classes" of quantum environment. 

We show how the mapping to spin bath models is made, and then discuss several examples in 
detail, including moving particles, magnetic solitons, nanomagnets, and SQUIDs, coupled to nuclear 
and paramagnetic spin environments. 

We then focus on the " Central Spin" model, which couples a central 2- level system to a background 
spin bath. It is the spin bath analogue of the famous "spin-boson" oscillator model, and describes, 
eg., the tunneling dynamics of nanoscopic and mesoscopic magnets and superconductors. We show 
how to average over (or "integrate out") spin bath modes, using an operator instanton technique, to 
find the Central spin dynamics. The formal manouevres involve 4 separate averages- each average 
corresponds physically to a different " decoherence" mechanism acting on the central spin dynamics. 
Each environmental spin has its own topological "spin phase", which by interacting with the phase 
of the central system, decoheres it- this can happen even without dissipation. We give analytic 
results for the central spin correlation functions, under various conditions. 

We then describe the application of this theory to magnetic and superconducting systems. Partic- 
ular attention is given to recent work on tunneling magnetic macromolecules, where the role of the 
nuclear spin bath in controlling the tunneling is very clear; we also discuss other magnetic systems 
in the quantum regime, and the influence of nuclear and paramagnetic spins on flux dynamics in 
SQUIDs. 

Finally, we discuss decoherence mechanisms and coherence experiments in superconductors and 
magnets. We show that a spin bath environment causes decoherence even in the T — » limit. Control 
of this decoherence will be essential in the effort to construct "qubits" for quantum computers. 



I. INTRODUCTION 



Many problems in quantum physics can be discussed using a model in which one or more mesoscopic or even 
macroscopic coordinates Q, Q' , etc., interact with a background environment (one cocujdinate might also represent an 
experimental probe, or even an observer). In such models (which have a long historyEI^H) all variables, including the 
environmental ones, are treated quantum-mechanically. The aim is to find the behaviour of Q, Q' , etc., after averaging 
over the environmental variables in some way. 

It is certainly not obvious that one can discuss the real world in this way, given the complexity of iV-body systems. 
However we now know that many (but not all) mesoscopic or macroscopic systems can be described at low energies 
by a few "canonical models", where a simple "central system" (eg., a 2-level system, or an oscillator) couples to an 
environment. Remarkably, there appear to be only two types of environment, describable as tjaths of either oscillators 
or spins. One way of trying to justify such models is the " renormalisation group" viewpointo, which maintains that 
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most physical systems fall into a few "universality classes", each scaling to its own "fixed point" in the space of 
possible Hamiltonians. All systems in a given university class use the same canonical Hamiltonian- the differences 
between different systems lying in the different values of the relevant couplings in this Hamiltonian. 

Although this "hard RG" philosophy clearly fails in some cases, it is a useful starting point for the present article, 
in which the quantum environment is modelled by a "spin bath" (usually of 2-level systems, or "spin-1/2" systems). 
The finite Hilbert space of each bath spin makes the spin bath appropriate for describing the Jf>w energy dynamics of 
a set of localised environmental modes. We concentrate on one particular "central spin" modeHu in which the central 
system itself reduces to a 2-level system; but we also discuss cases where the central system is a one-dimensional 
coordinate (a "particle") moving through a field of spins. 

Another well-known set of caafitfiical models describes the envirompent as a set of uncoupled oscillators- these 
include the " spin- boson" modeOO and the " Caldeira-Leggett" modcO'E^I. The spin-boson model couples a central 
2-level system to the oscillators, and is thus the analogue of the central spin model; and the Caldeira-Leggett model 
couples a timpeling particle to the oscillators. These oscillator models all derive from a scheme proposed by Feynman 
and Vernont3, to describe a central system coupled weakly to N environmental modes; as they showed, the mapping 
to an oscillator bath can only be made rigourously if the coupling is weak. Oscillator models are thus best adapted 
to N delocalised environmental modes (where the coupling is automatically ~ l/A^^/^, and thus small for large N). 

However, readers familiar with low-temperature physics will know that at low energies, the entropy and heat 
capacity of almost all real physical systems are dominated instead by local modes such as defects, impurity spins, and 
nuclear spinscj. Typically these relax very slowly at low T because little phase space is available in their coupling 
to any delocalised modes (or to each other). However they often couple strongly to any mesoscopic or macroscopic 
collective coordinate, which then easily perturbs them. This coupling is of course independeni. of N. 

Unfortunately, even though spin bath models have been studied sporadically for marui-iitoarsS, the results have often 
been misleading, either because they treat some weak-coupling limit (sometimes madet^lla by arbitrarily multiplying 
the coupling to each of the N bath spins by a factor for no good physical reason), or because they drop some of 
the important couplings to the batJji|Si)ins, in order to solve the model. In the weak-coupling limit, spinJaath models 
can be mapped to oscillator bathalZl^lij (in accordance with the original remarks of Feynman and VernonEJ). 
one is often nowhere near the weak-coupling limit, and the mapping to the oscillator bath then fails in genera. 
This demands a new approach, which is the subject of this review. 

It may be useful to mention why many physicists are interested in models of this kind. Here are some of the reasons: 

(i) Very rapid progress in work-aiLintrinsically quantum processes (interference, tunneling, etc.) occurring at the 
nanoscopic and mesoscopic scalesEEl cJ, plus speculations about the coming " nanotechnological revolution" . Perhaps 
the most excitiag idea in this area is that of making "quantum computers" using nanoscopic superconductorsEj, 
semiconductorstZl, or nanomagnets. Needless to say, the technological repercussions of this work will be enormous, 
provided the crucial problem of decoherence can be overcome. ■— ■ 

(ii) Physicists need to understand the mechanisms of decoherence and quantum dissipation^ in nature, and the 
crossover to (or "emergence" of) classical behaviour from quantum physics as either size, temperature T, external 
fields, or couplings to the environment are increased. These issues are not only relevant to quantum device design, 
but also to problems in quantum gravity, and to the infamous "quantum measurement" probleinc3. The existence of 
low-T canonical models, going beyond the phenomenology of stochastic or master equationtO to work-with closed 
Hamiltonians, is invaluabla-bjO-L-p ..Rccf-nt examples include the analysis of uuHntyLm-spin glass relaxations, quantum 
relaxation in nanomagnetsBu ljEjEj'LiJ- iJ and tunneling superconductorst3E3 ii cJoo) , and the study of quantum 
chaotic systems, stochastic resonanceea and dissipative tunnelingJjji AC fieldsCJ. Earlier such models have been used 
for decades to discuss relaxation in fields like quantum chemistr5i£3 or nuclear physics. 

(iii) Both oscillator and spin bathj-inodcls map to many important models in quantum field theory. Thus the 
"spin-boson" model mentioned abovellB maps, for specific parameter values, to the Kondo model, the Thirring and 
Sine-Gordon models, and various other 2-dimensional field theories. Although similar mappings have yet to be 
exploited in great detail for the spin bath, they will obviously be very useful for, eg., lattice spin models. 

Most work in these areas has used oscillator bath representations of the environment, with the tacit assump- 
tion thai—delapalised environmental modes dominate the physics. However, experiments on quantum nanomagnetic 
systemsOEJcJcd, on glasseaEl, and on mesoscopic conductorsEjS, as well as theoretical debate about the mechanisms 
of decoherence in nature, clearly require a more general point of view. In fact we shall see that spin baths behave 
very differently from oscillator baths. For example, the oft-discussed connection between dissipation and decoherence 
which, exists for oscillator bath environmentsc3'Eil is absent here- one can even have decoherence with no dissipation 
at allEI, because of the quantum phase associated with the spin bath dynamics. 

Although the existence a£.a quantum phase associated with spin is obvious (it is a quantum variable), it was not 
until HaldaneL3 and BerryE3 discussed its topological properties that physicists realised its practical importance, in 
terms of the path traced out on the spin sphere. Just as in the usual Aharonov-Bohm effect, the "flux" enclosed by 
a path (given here by ujS, where S is the spin and ui is the enclosed solid angle on the unit spin sphere) is equal to a 
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dynamical phase- but now these are both in spin space, not in real space. These ideas (and related experiments) were 
the eeptre of enormous interest in the mid-late 1980's, in almost all fields of physics (and were extensively reviewed 
thenN). 

In this article we will be interested in the spin phase of the environment. We stress here that the environmental 
"spin" variables may not necessarily refer to real spins (they can describe defects, or other such "2 level systems"), 
but they will still have an associated dynamic topological phase, which can be described by an effective spin bath 
variable. The.|Cnjfironmental spin phase interacts with the phase of the central system, causing phase decoherence in 
its dynamicsElBO. From the point of view of measurement theory, this environmental phase decoherence comes from 
a "phase measurement" made by the spin environmenttll, in a kind of "inverse Stern-Gerlach" setup (where the spins, 
instead of being measured, are themselves doing the measuring!). Such phase decoherence also arises from oscillator 
baths, in a rather different (and much less effective) wayE3. In fact, the relevant phases involve both an adiabatic 
"Berry" term and a second term coming from transitions between different bath states (section 3. A). There are also 
other decoherence mechanisms associated with the bath spins, coming from both the temporal fluctuations of the bias 
on the central system caused by the spin bath (section 3.D), and from the precession of the spins in the spin bath 
(with their associated phase) in between transitions of the central system (section 3.C). Thus the question of how the 
bath spin dynamics influences the central system is not simply a question of looking at Berry phases. 

Practical application of the theory to, eg., SQUIDs, or nanomagnets, or "qubits" (section 5), must include |-aU. 
mechanisms properly (section 3.E). The tactic adopted in this article is to focus on a "Central Spin" modeMtl 
(sections 3 and 4), in which the role of each term is exposed rather clearly; after this one sees how things generalise 
to other models. This model is directly relevant to qubits, and to the observation of mesoscopic or "macroscopic" 
quantum coherence- indeed we maintain that any practical design of such devices, must involve the elimination, by 
one means or another, of all decoherence mechanisms from the relevant spin bathtJ. 

We begin the article (section 2) by showing how both oscillator and spin bath models arise as the low-energy 
truncated versions of higher energy Hamiltonians. We give several examples, both magnetic and superconducting, to 
illustrate this. Then, in section 3 we explain how, mathematically, one averages over the spin bath variables to find 
the behaviour of the central system. This is done pedagogically- we use the example of the Central Spin model (and 
compare it with the spin-boson model). Various simple limits are introduced, and solved, before the general technique 
is given at the end of section 3. In Section 4 we give some resulta] for the dynamics of the Central Spin model, in 
various regimes, to show how the physics is influenced by the spin bath; and we also show how the system reduces to 
the spin-boson system in the weak-coupling limit. For those readers interested in the mathematical details, these are 
sketched in 2 Appendices. 

Finally, in section 5 we return to physical applications, particularly to quantum magnetic systems and superconduc- 
tors. We then discuss decoherence, and show how this should persist even in the T ^ limit. We finish by discussing 
the important application of the Central Spin model to decoherence in qubits and in quantum computation. 

II. THE LOW-ENERGY EFFECTIVE HAMILTONIAN 

In studying the low-energy dynamics of a central quantum system coupled to an environment, we begin by "trun- 
cating out" the unwanted high-energy physics, to produce a low-energy effective Hamiltonian. This is of course a 
quite general technique in physics, and one way to approach it is illustrated in Fig.l. Typically one has a reasonably 
accurately known "high-energy" or "bare" Hamiltonian (or Lagrangian) for a quantum system, valid below some 
"ultraviolet" upper energy cut-off energy Ec, and having the form 

HB.rc{Ec) - Ho{P,Q) + Hint{P,Q;p,q) + Henv{p,q) {E < E,) , (2.1) 

where Q is an M-dimensional coordinate describing that part of the system we are interested in (with P the cor- 
responding conjugate momentum), and {p,q) are A^-dimensional coordinates describing all other degrees of freedom 
which may couple to {P,Q). Conventionally one refers to {p,q) as environmental coordinates. i?Barc is of course 
a low-energy form of some other even higher-energy Hamiltonian, in a chain extending ultimately back to quarks, 
leptons, and perhaps strings. 

If, however, one is only interested in physics below a much lower energy scale ilo, then the question is - can we find 
a new effective Hamiltonian, of form 

H^sirio) = Ho{P,Q)+Hi,,tiP,Q;p,q)+H,,-UP,q) {E < n^) , (2.2) 

in the truncated Hilbert space of energies below SIq? In this Hcs, P and Q are generalised m-dimensional coordinates 
of interest, and p, q are A^-dimensional environmental coordinates coupled to them. Since we have truncated the total 
Hilbert space, we have in general that M < M and N < N. 
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Why do we make this truncation (after all, its inevitable effect will be to generate new couplings between the 
low-energy modes)? Essentially because in many cases the truncation pushes the new H^s towards some low-energy 
"fixed point" Hamiltonian; and many different physical systems may flow to the same fixed point. This allows us to 
speak of "universality classes" of quantum environment, and of a small number of "canonical" effective Hamiltonians. 
All physical systems in the same universality class will be described by the same form for i?cff, albeit with different 
values for the couplings. As one varies the UV cut-off Qo, the couplings change and any given system moves in the 
"coupling" or "effective Hamiltonian" space; but they all move towards the same fixed point (or fixed line, in some 
cases) as fio is reduced. The various coupling terms in H^f f{^o), simply parametrise the path it takes as llo — > 

Explicit derivations of i/cff for particular systems are leBgt]w;|-se&|es£llll'il'ii'ii'S for general discussions, ref.El 
for comparision of spin and oscillator bath systems, andt3t30'cZlBl23'Ej for specific examples. In this article we 
will go directly to the canonical models, giving some examples of each so that readers can see how the high-energy 
Hamiltonians are related to the models for some real systems. To warm up we recall the basic structure of the 
oscillator bath effective Hamiltonians, and then move on immediately to discuss various canonical models involving 
spin baths. 



A. Oscillator Bath models 



For models in the general "universality class" of oscillator bath environments, iJcff takes the for 



nJiS: 



N 



Fg{P,Q)Xg+G,{P,Q)p, 




Pq , 2 2 

1- niqUJqXg 



(2.3) 



where the generalised bath coordinates {qk,Pk) are now oscillator displacements Xq and momenta p^; these describe 
delocalised modes. The couplings Fq{P,Q) and Gq{P,Q) are ^ 0{N~^/'^), so that in the " t herrn nd 3m amic limit" 
^ 1, appropriate to a macroscopic environment of delocalised oscillators, these couplings are smalLj'E^EJ (a nptaber 
of studies have also shown how higher-order couplings can be absorbed into linear but T-dependent couplingst2l) . A 
special case of (^.3[) is the Feynman- Vernon bilinear coupling formlij: 



H,aino) = Ho{P,Q) 



N 
9=1 



CqXqQ 



TUqUjlxl 



(2.4) 



9=1 



where the couplings Cq ^ O^N" ' ) as well. In recent years great jatteiition has been given to problems where 
Ho{P,Q) describes a tunneling system (the " Caldeira-Leggett" modeOO); thfiC^ have also been extensive studies 
of an oscilLatiir coupled to oscillatora£J, of free particles coupled to oscillatorS^H and of band particles coupled to 
oscillatorsU'El. 

Suppose now the potential V{Q) has a 2- well form, with small oscillation frequencies ~ flo, and a "bias" energy 
difference between the two minima < flo- Then for energies < flo, one further truncates to the celebrated "spin-boson" 
modetj: 



N 

E 

9=1 



h.C.) 




Pq , 22 

h niqUJgXq 

9 ; 



(2.5) 



where the two-level central system (with tunneling amplitude A($o) and longitudinal bias ^h) is described by the 
Pauli matrix vector t, coupled to background oscillators having energies LUq < Qq- We have introduced a topological 
phase $0 for the central system, which depends in general on an external field; the simplest and best-known example 
is a form A(<I>o) = 2AoCOS<I>o, arising from the interference between 2 paths of amplitude AqC***" in the motion of 
the central system. This kind of " Aharonov-Bohm" interference is well-known in superconductors (where theLflhase is 
just the superconducting order parameter phase), and in magnets (where it is the topological spin phaseE^TEJ). One 
can have a more complicated dependence of A($o) on $0 (eg., using multiple SQUIDs, or a nanomagnet with ni-fold 
rotation symmetry), but we will stick with the simple 2cos$o dependence in this article. 

For consistency we must assume £,h < ^o, otherwise higher levels will be involved. Typically is dropped, because 

its effects are down on those of c| by a factor ^ (Ao/ilo)^ in tunneling rates; but sometimes cf = (for reasons of 



symmetry) , and then must be retained. The tunneling amplitude Aq ^ 



where is the tunneling action. 
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The fame of the spin-boson model partly arises because many well-known problems in condensed matter physics can 
be mapped to it- this is a good example of the "universality" mentioned above. Because the effect of each oscillator on 
the central system (and vice- versa) is very small, it may be e ntirely incorporated in second-order perturbat ion theory 
(ie., to order ~ (F^/ujq), (G^/w,)) for the general form (2^), or to order ~ {c^/ujq) for the bilinear forms (2.4),( ^!5| ). 
In the latter case this immediately encapsulates all environmental effects in the spectral functionHEj 



Jaiuj) 



N 

2 ^ 

9=1 



(2.6) 



bqujq 



where a =_L, ||. In general Jai(^) also depends on T, even in the \ow-lu limitEiltfi'3. The case where Jq(w) ~ is 
referred to as "Ohmic"ll3. Because the Cq ^ N^-^^^, the Jq(w) are independent of N and have the usual "response 
function" form. 



B. Two examples of Spin-Boson systems 

We give just 2 examples here of how a spin-boson model can arise, in the description of mesoscopic systems at 
low energies. Both truncations ignore the presence of spin bath modes (for which see sections 2.E and 2.F, where we 
return to these 2 examples). 

(i) Nanomagnet coupled to phonons or electrons: The electronic spin dynamics of nanomagnets are often 
described by-a^giant spin" Hamiltonian Ho{S), describing a quantum rotator with spin quantum number S = \S\ 3> 1. 
This modelEjO assumes the individual electronic moments are locked together by strong exchange interactions Jij 
into a monodomain giant spin, with S = sj (summed over local moment sites) . This only only works0 below a 
UV cut-off energy Ec considerably less than Jij . However we are interested in the quantum dynamics for energies 
< ilo, where flo is controlled mostly by the single- ion magnetic anisotropy; in real nanomagnets Jlo ~ 0.1 — 10 K. 

Any real nanomagnet has couplngs to a spin bath of nuclear and paramagnetic spins, and to oscillator baths of 
phonons and electrons, which we now describe. i— ii— i 

The " high-energy" coupling between phonons and S is described by terms likeE3E3: 

where mg is the electron mass. Ma the mass of the molecule, and ^ c^a"^ is the Debye temperature (with a 
the relevant lattice spacing, and Cs the sound velocity). The interaction U{S) ~ S and dimensionless; a typical 
example is the non-diagonal term (SxSz/ S), which causes phonon emission or absorption with a change ±1 in Sz 
(since Sx = ^{S+ + S^)). One also has diagonal terms in which SzSx is replaced by, eg. S'J; and there are also higher 
couplings to, eg, pairs of phonons. j-. 



Truncation to the "quantum regime"l-l then gives the spin-boson model (2.5), with a dominant non-diagonal coupling 
^ STJolffl^^^, coming from terms like (p/7|). In the absence of external fields in Hamiltonian H^^-f (S), the diagonal 

coupling Cq is actually zero (because of time-reversal symmetry). The Caldeira-Leggett spectral function for the 
system has the form J±{uj) ~ {^l^ / pc^)uj^ where B± ^ (S'^il^/Oi)); here p is the density of the medium supporting 
Debye phonons, and 8|) ^ pc^. P 

If now we couple electrons to the giant spin, it is the diagonal coupling which dominatesEI. The electronic coupling 
to S depends on the type of magnetism. Some details have been worked out for Kondo interactions with conduction 
electrons- the coupling to S is 

kq 

where J is the mean value of the Kondo couplings to each individual electronic spin in the nanomagnet, and SFq = 
J {(Pr /Vo)s{r)e'^'''^ is a "form factor" integrating the localised electron spin density s{r) over the nanomagnetic volume 
Vo- At low energies the corresponding spin-boson model has an "Ohmic" diagonal spectral function J\\{uj) — Tra^^. 
The size of depends on how the electrons permeate the nanomagneta; if they permeate freely. 
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where g = JN{Q) is the mean dimensionless Kondo coupHng, and A^(0) the Fermi surface density of states. Typically 
g 0.1, so Q!k can be big. 



(ii) RF SQUID (Flux coupled to electrons): We briefly recall one well-known application of the si 
model, to an RF SQUID coupled to both normal electrons (in shunts, etc.), and Bogoliubov quasiparticlca 
The flux (j) passing through a superconducting ring with a weak link moves in a multiwell potential, which can be ad- 
justed so that the lowest 2 wells (each with small oscillation or " Josephson plasma" frequency Qo ~ '2tt[Ej /nC]^^^ /(j)o, 
where Ej is the Josephson weak link coupling energy) are almost degenerate, and dominat|e-.thft low-energy properties. 



The high-energy coupling between the flux and the electronic quasiparticles has the forr 

H^nt = {cos(7r<?!>/0o) UUq{aq + 4) + i sin(7r0/,^o) ^ tqUf^a^ - aj)} (2.9) 

9 9 



where q = {k,k') labels oscillator states describing a quasiparticle pair \kk') with energy ujq = + E'^/, tq is the 

relevant junction tunneling matrix element, Uq^^ the symmetric/antisymmetric BCS coherence factor, and (j)o is the 
flux quantum. Thus we have a coupling to both the morpenta and coordinates of the oscillators, which can also 
be written as a coupling to 2 independent oscillator bathscJ. The T-dependence of the coherence factors (coming 
from the BCS gap dependence) as well as the gap structure in their energy dependence, gives a complex structure 
in J{llj,T). The reduction to the spin-boson model is now triviaO, the minima in 0-space of the effective potential 
corresponding to the 2 eigenstates of . 



C. Spin Bath Environments 



Now suppose we have a high-energy Hamiltonian of form (2.1), but where the environmental coordinates {p,q) are 
a set of N spin - 1/2 variables {f?fc}, (i.e., two-level systems); and we assume the interspin couplings to be weak. Then, 
instead of (2.3), we have 



H = Ho{P, Q) + -ffi„t(P, Q; {S}) -f i?env({5=}) 



(2.10) 



N r 



i?int(P,Q;{c?}) 



fc=i 



(2.11) 



N 



N N 



k=l 



k=l k' = l 



(2.12) 



for energy scales E < Ec- Thus we now have a central system coupled to a "spin bath", described by Hcnv{{<^}) 
in ( ^.12 ). The couplings f|(P, Q) and Fj^{P,Q), between the central system and the bath spins, are usually much 
greater than the interspin couplings V^j^ ; this means that the dynamics of each spin is largely " slaved" to that of the 
central system. 

Unlike oscillator baths (whose modes typically represent delocalised environmental degrees of freedom), the {(Jk} 
represent localised modes (whose weak spatial overlap explains why the V^/.^ are small). This fact underlies a crucial 

difference between oscillator and spin bath environments- the couplings i^|(P, Q) and Fj^{P,Q) are independent of 
the number N of bath spins. Thus the larger is TV, the larger is the total effect of the spin bath on the central system- 
there is no strict thermodynamic limit in the system, and it is not meaningful to let ^ oo. We emphasize also that 
we see no justification in general for spin bath models in which F^^jFj}- '-^ 0(A^~^/^), or even ~ 0{1/N) (although 
one can certainly invent artificial models of this kind). Thus, if we add more localised environmental modes to our 
environment, it is clear that the different modes are approximately independent (as they will be if quasi- localised), so 
that their individual couplings to the central system will be hardly affected, ie., will depend only weakly on N. 

There is nothing to stop generalisation of this model to include bath spins {Ik}, with 1^ = \Ik\ > 1/2; the (2/^ -I- 1) 
states then represent the degrees of freedom. of, eg., a defect, or a spin (again, localised). This introduces tensor 
(eg., quadrupolar) couplings to the bath spinsEjlZ3, and thereby complicates the mathematics- but does not alter the 
basic physics. We wU]-|Bot discuss this here (for the relevant formalism, and its application to the Fe-8 molecular 
nanomagnet, see refsHj'E2l). 
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D. Particle moving through a spin bath 



A particle moving through a spin bath is described by ( 2.10| ), in which P and Q describe the momentum and 

II 

position of the moving particle. The diagonal term Fj^{P,Q) is analogous to the "position" oscillator coupling 



Fq{P, Q) in (2.2), and likewise F^{P, Q) to corresponds to Gq{P, Q). However both forms can be altered by canonical 
transformation, corresponding to a rotation between the different coordinates. The most common problems involve 
a diagonal coordinate coupling -F'I(Q) and a non-diagonal momentum coupling F^{P). Then bath transitions (spin 
flips) are induced by the motion of the particle, whereas a stationary particle sees a "potential" [/(Q, {(T|(t)}) = 

Fj,'(Q)(T|(t), in general time-dependent. r-, 
A nice mesoscopic example of this is a large magnetic soliton coupled to background spinsE3. In many realistic 
cases the most important such coupling will be to paramagnetic impurities, but here we consider the simpler case of a 
hyperfine coupling to a set {t?} of N spin-1/2 nuclear spins. In this case Vl^^ describes the extremely weak internuclear 
dipolar couphng; typically \Vl^^, \ < 10^^ K; and hk is any external field that might unfluence these nuclei. 

The "high-energy" HamiltonianJar such a wall is usually determined as an integral over the magnetisation density 
M(r) and its spacetime gradientalJ. ^From this one eliminates the details of the wall profile altogether, to pro- 
duce a "bare" Hamiltonian (ie., neglecting the environme nt j) fo r the wall coordinate; in simple cases where the wall 
demagnetisation field keeps the wall almost flat, this givescilO 



-M^Q^ - V{Q) - 2S^iibMoH,Q 



(2.13) 



for a wall with centre of mass coordinate Q and surface area Syj. The "pinning potential" V{Q) — Vq sech^((5/A^), 
provided the pinning centre is much smaller than the wall width A^. The term linear in Q comes from an external 
magnetic field Hg. 

WhaL.p£iw of the environment? IrLlhe literature there is extensive discussion of the effects of magnons (ie., spin 
waves)EjO, electronsO, and phononsEJ on the wall dynamics- these are all oscillator baths. However at low T spin 
bath effects, coming from nuclear and paramagnetic spins, will completely dominate. Even in Ni (where only 1% of 
the nuclei have spins, with a tiny hyperfine coupling = 28.35 MHz ~ 1.4 mK), all real samples have an important 
concentration of paramagnetic spins (caused by Oxygen in the sample) as well as many defects. In rare earths, the 
hyperfine coupling LOk ^ 1 — 10 GHz (0.05 — 0.5 K), and hyperfine effects alone are quite massive. Thus we must 
modify above to 



N 
fe=l 



a/3 „a tP 



Q/3 TO. Tp 

kk' ^k ^k' 



(2.14) 



in which the electronic spins couple locally to N nuclear spins 1^: at positions {k = 1, 2, 3, ... N), via a hyperfine 
coupling w^'' (and also in general to paramagnetic spins). The internuclear coupling \Vl^f^\ ^ 1 — 100 kHz (0.05 — 5/iK), 
ie., <ti LUk, but it gives the spin bath its o wn d ynamics. 

To write the Hamiltonian in the form (2.11), we write the continuum magnetisation M(r) = Mo(r) + m(r), where 
Mo(r) is the slowly-varying part descibing the wall profile and m(r) clescribes fluctuations around this. Then we 
rotate the spin quantisation axis to be locally parallel to Mo(r), and getE3 



N „ ,3 r 

H = H.^+Y, — <5(r-r,) c.j!Af,(r)4^ + ui[m^{v)Il + my{v)Il] 
k=i-^ ^9 L 



^EE 

k k' 



y^k^ikit^ (2.15) 



displaying explicitly the longitudinal and transverse couplings. The "particle" moves through a slowly fluctuating 
"random walk" potential field t/||(Q) (coming from the sum over couplings wj! to randomly oriented spins). The 
transverse coupling (independent of Q but not of P) causes "spin flip" transitions in the spin bath when the wall 
moves, even if the bath is at T = 0. 

One may also discuss problems in superconductors and normal metals involving nuclear and paramagnetic spins, 
and other " defects" , which can also be mapped to the same model of a particle moving through a spin bath (sections 
2.F, 5.B). 
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E. The Central Spin Model 



Is there a low-energy effective Hamiltonian, analogous to the spin-boson model, in which a " central" 2-level system 
couples instead to a spin bath? The answer is yes, but the effective Hamiltonian does not look qui te so simple as the 
spin-bosoH iQp^. In the absence of any external field, the analogue of the spin-boson form in (2.5) for a spin bath is 
actualh 



2Af_ cos 



N 

fc=i 



N 



H.c. 



N N 



(2.16) 



k=l 



k=l k' 



where r describes the central spin, and the Uk the spin bath degrees of freedom. This form is not the most general 
one- apart from dropping external field effects (for which see below) we have also restricted the central spin phase 
in a simple cos$ form (c/. Introduction). As discussed below, both A and $ incorporate spin bath renormalisation 
effects. The factor of 2 in front of A is somewhat arbitrary (if the cosine is one, then the actual "tunnel splitting" 
coming from ( ^.16 ), i n the absence of spin bath effects, will be 4A). 

The basic form of ( 2.16| ) is actually fairly easy to understand. The extra phase in theJtirst "non-diagonal" term 
(adding to $) comes from the topological phase of the bath spins as they make transitionsH. There are also diagonal 

terms, and a weak interaction Vkk' between the bath spins. This form assumes the diagonal couplings a;|,a;^ are <C 

the UV cutoff Q.O, and that Vkk' <C uj\,uj-^. The ratio Vkk' is arbitrary. In sections 4 and 5 we shall see there is a 
weak coupling limit to this model, in which it reduces to a spin-boson system. 

Before this model was derived in this general form and then solved, a number of special cases had already been 
looked ati. In particular, Shimshoni and Gefer0 included only|ptiw_ diagonal terms aj| and w^, and examined the 
results in weak coupling in the presence of an AC field (see alsoESCj'Ej) ; they clearly recognised that the problem was 
different from an oscillator bath one. 

Let us now discuss the different terms in (2.16), in the order they appear (cf. Fig. 2). 



(i) Non-Diagonal terms: That an extra phase term should exist, coming from the spin bath, is obvious on general 
grounds (cf. Introduction). One can understand its algebraic form in the following way. Notice that -ffe// in ( 2.1(: ) 
operates on both the central spin and the spin bath; and the effect of a single giant spin transition on Sk can always 
be written as a transition between initial and final nuclear states, in the form | = Tk \ X™), where 



. — e J 



(2.17) 



The integral over Hint is onlyjiefined once we know the trajectory of the central system during the tunneling event. 
It is in this sense that we sayO that the instanton has become an "operator" in the space of the spiruhalh modes. 
Notice that in general the central spin phase $ and splitting A are renormalised by the bath couplingsEiL 



^ = ^o + ^<i>k 



A = Aoexp{-^(5fe} 



(2.18) 



The {(pk} are "Berry phase" terms coming from the bath spin dynamics during a central spin flip, and the {5k} come 
from high-frequency modifications of the original high-energy potential ("barrier fluctuations"). One can show that 
both (j)k and 5k are ~ 0{uj'^/Qj'^), where ujk is the larger of and we will ignore these terms from now on. 

Expanding out the cosine in ( ^.16| ) gives a series of terms like 'T'±^ ap-y---(^k^c^k2^'k3 • • • in which the instanton flip 
of the central spin couples simultaneously to many different bath spins- a single central spin transition can cause 
multiple transitions in the bath (Fig. 2). Later we will introduce a parameter A which measures the average number 
of bath spins flipping during each instanton. 

(ii) Diagonal terms: These act between transitions of the central spin (Fig. 2), and are also easy to understand. 
Formally, one starts by considering the "initial" and "final" fields (ie., before and after a transition of r) acting on 
fjk- Calling these fields 'f^^^' and 7^^'' respectively (Fig. 3), we define the sum and the difference terms as 



(1) 



t^fe"m, = (7r+ 7^)72. 



(2.19) 
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where Ik and rhk are unit vectors. Then the truncated diagonal interaction takes the form 



= Y.{€'-^+l'r-^\-^>^ ^ r.J24h-lfk+j:^irnk-h, (2.20) 

fc=l ^ ^ k=l k=l 

i.e., one term which changes during a transition of the central system, and one which does not. 

The longitudinal coupling Tz'^Luj.a^ determines the gross structure of the bath states in energy space: it also 
determines an "internal bias field" e({cr|}) — J^^^k acting on t^. The 2 levels of bath spin dl are split by energy 
wj!, depending on whether cr| is parallel or antiparallel to t^. The effect of this on the 2^ fold multiplet of bath 
states surrounding each central spin state is shown in Fig. 4. Suppose we classify these states by their "polarisation 
group"; all bath states whose total longitudinal polarisation J2'^k — i^- polarisation group M. Since the cj|. 

vary from one bath spin to another, states in polarisation group M are spread over an energy range Tm ] and the 
entire manifold of states, comprising all polarisation groups, is spread over a larger energy range Eq- Let us define 
normalised densities of states Gi\i{e) and W{e) for these 2 distributions, so that 



Wie) = (1/2^) J2 c(,™^/'GM(e) (2.21) 



where C™ = nl/ml{n — my.. In almost any physical case one will have strongly overlapping polarisation groups, so 
that for all but very small values of N, or except in the extreme wings of the distributions, one has 

GM{e) - (2/7rf2,)V2e-2eVfL_ (2.22) 

T^(e) ^ (2/7ril;2)l/2e-2^Vi^„^ (2.23) 
The simplest case is where the wji cluster around a single central value lUq with variance 

= ^I^TM^' (2.24) 

For this case we define a parameter /i = 7V^/^5wo/wo, characterising the degree of polarisation group overlap; overlap 
is complete when /z > 1. Then Tm ~ 2N^/'^SuJo and Eo = 2N^/^uJo (so that Tm/Eo ~ Suo/uio = N'^^'^^i). In 
the extremely unlikely case where /i <C 1, W{e) can no longer be treated as Gaussian- however there is an intrinsic 
lower limit to the linewidth of each polarisation group set by the interspin interaction Vkk' ■ This " intrinsic linewidth" 
To ~ iVi/^K, where Vo is a typical value of Vkk' ; in any physically realistic case this is usually enough by itself to 
cause complete overlap of all groups (the essentially non-interacting case where F = 0, ie., ujk — lOq for all bath spins. 



and Vkk' = 0, 14 = 0, lo^ = 0, was actually studied by GargEll 

The "transverse" couplings Lu-j: arise when the fields before/ after a transition, acting on the {f?fc}, are not exactly 
parallel or antiparallel. This can happen in many ways, either because of external fields which couple to the bath 
spins, or because of a lack of symmetry in the underlying dynamics of the central system, or in its coupling to the 
bath modes. Thus they are non-zero in any realistic situation. 

(iii) Internal Spin Bath dynamics: Finally, the interaction Vkk' is usually so weak that it does not change under 
truncation. If this term is absent, the spin bath will have no "intrinsic" dynamics, and remains inert between 
transitions of f. Thus even if small, Vkk' is important, since it allows the bath state to evolve during these intervals. 
The most important effect of Vkk' is that it allows the longitudinal bias field e({cr|}) tofluctuate in time, between 



and during transitions of r. Notice, however, that with the Hamiltonian in the form (2.16), only fluctuations within 
the same polarisation group are allowed. In NMR language, only T2 processes occur in the intrinsic dynamics of the 
spin bath- changes in M can only occur via the interaction with the central spin. This will only be true at low T- at 
higher T longitudinal relaxation (ie., "Ti processes", in NMR language) between different polarisation groups should 
be included in a realistic model. Such processes (which arise from the interaction of the spin bath modes with other 
environmental modes, or with thermally excited higher modes of the central system, above the UV cutoff energy J7o) 
are almost always very slow when kT <C fio- On the other hand T2 fiuctuations will persist until kT ^ Fq; in the 
physical examples studied so far this means they persist down to temperatures or below. We will return briefiy 
to this very low T regime at the end of the article (section 5.C). 
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External Field Effects: In general all of the parameters in (2. If) will depend on any external field Ho, because 



it changes the high-energy dynamics of both central system and bath- however we can make low-field expansions and 
separate out the most important terms. Defining the "Zeeman" coupling energies fln^ and u^" of central and bath 
spins to this field, it is easy to see that under the conditions fln^/i^o < 1 and uj^° /ujk < 1, the principal changes to 



( p. 16 ) will be (i) the addition of an obvious longitudinal coupling S.Ho^'z to the central spin, and (ii) Qthe changes 
$o^*ff„, Vfc->y;-"°,and ujj: = Lo^rhk a;^(Ho), where up to linear order in Ho one has 

= 14 + vu{Ho) ; \vk{Ho)\/\Vk\ - 

L^tiUo) = Lji + dfe(H<,) ; \dk{llo)\/uji ^ c^f°M- (2.25) 

Thus even at low fields one has an important change in all the topological phases in the problem, and also to the 
transverse diagonal coupling (which itself arises from internal fields). In general Hq will also change the interspin 
bath couplings Vkk' in a way which depends on the specific details of the problem. In the next 2 sections we see how 
this works for both magnetic and superconducting systems. 



F. Nanomagnet coupled to nuclear and paramagnetic spins 

If we start from the "giant spin" model introduced above for a nanomagnet, then a simple isotropic contact hyperfine 
coupling to nuclear spins will lead to a Hamiltonian (for E < Ec) like: 



H{S;{<j})^HoiS) 



1 



N 



^ k=i 



(2.26) 



where Ho{S) is the "giant spin" Hamiltonian, and i?onv({(?}) is the same as in (III). The generalisation of this 
simple Hamiltonian to include dipolar hyperfine interactions, as well q^^Q higher spin nuclei and to paramagnetic 
spins (with tensor and quadrupolar couplings) can be used if necessaryEj£3. Here we will assume for simplicity that 
I Ik 1= / = ^, and write Ik cfk, i.e. the nuclear spins will be described by spin-i Pauli matrices. In fact in many 
cases even ii I ^ ^ , the low-energy nuclear spin dynamics is well described by a 2-level system. 

The ni^ff¥ p?f^^f*^P '^^ ^i^ii'^}) to a central spin Hamiltonian 7foff(r; {c?}) has been discussed in several 
papersBLta'Prr^r i. As an example we quote the result for a simple easy axis-easy plane nanomagnet (for which 

HoiS) — {1/S)[~kIs'^ + K2Sy], and give it a phy^cj^^ij^nterpretation. In this case, assuming tok <C and also a 

weak external field Hq, the effective Hamiltonian is 



Heffi^o) = < 2Aof_ COS 



vrS* - i ^ akfik ■ Sk - /^oHo-Hc 



N 



fe=l 



N N 



k=l k' = l 



H.c. 



(2.27) 



ie., a special case of the general form (2.16), with the parameters = gt^-BSzH^, h = 



= u!k, and uj-ir = 0. The 



vectors a^rifc and (3o^o for this easy axis-easy plane case turn out to be (again assuming ujk ^ 51oj and small Hq): 



aknk 



(2.28) 



In this example there are 2 tunneling trajectories (clockwise and counterclockwise in the easy plane), giving a result 
giiQfcrifc o-fc ^Yic "transfer matrix" Tk (in zero applied field). The resulting vector akfik is the "average hyperfine 
field" acting on cjk during the tunneling event. To understand its orientation (and why it is complex) we note that 
the nuclear spin itself exercises a torque on S while it is tunneling, and this pushes S away from the easy plane. 
Consequently (a) the average field acting on Uk has a component out of the easy plane, in the y-direction, and (b) 
S no longer moves exactly along the easy-plane path, while tunneling, that it would in the absence of dk (and so its 
action increases, via the imaginary part of akfik)- 
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The contribution /^qIIo.Ho to the topological phaseEj'LatZl comes from the area swept out by the giant spin on 
the spin sphere (cf. Introduction), which changes, as the field Ho changes; it is essentially an " Aharonov-Bohm" 
contribution to this phase from the external fieldO, which leads to spectaculax_Qacillations (Fig. 5) in the effective 
tunneling amplitude A = 2Ao cos[7r5' + i/3ono.Ho] for a field perpendicular to zOEZl. Very reeently oscillations in the 
tunneling amplitude of Fe-8 magnetic molecular crystals were seen which are related to thitfia, although the presence 
of both nuclear spins and dipolar fields seriously complicates their interpretatior£3 (see section 5. A). 

We recall from section II. B that the giant spin model truncates for nanomagnets to a 2-level system for energies 
<C r^o ~ 1 — 10 if . Contact hyperfine couplings are in the range 1.3 — 30 mK (transition metals) or 40 — 500 mK 
(rare earths); on the other hand the internuclear couplings Vkk' ^ 10^^ — 10^^ K. In the Fe-8 system just mentioned 
the hyperfine interactions are actually dominated by dipolar couplings between the 8 Fe^^ ions and the 120 protons 
in the molecule; these couplings are in the range ^ 1 — 100 MHz (0.05 — 5 mK). When the hvnerfine couplings are 
this weak we must also take into account the effect of external fields on the nuclear dynamicalJ (section 5. A). The 
values of A vary over a huge range, but typically A ^ ujk (in Fe-8, A ^ 10~^ K). 



G. SQUID coupled to nuclear and paramagnetic spins 

We consider again the RF SQUID, but now concentrate on the coupling of the flux $ to the spin bath of nuclear 
and paramagnetic spins which are within a penetration depth of the surface of the superconductor. This example is 
very instrpftive in understanding the weak-coupling limit of the central spin model (the following discussion is based 
on refs.NO). j-. 

Suppose to start with we consider a "cubic geometry" ta, in which a cube of side L — 1 cm has a hole of radius 
R — 0.2 cm through it, with in addition a slit connecting the hole to the exterior, spanned by a cylindrical junction 
of length I — 10~^ cm and diameter d = 2 x 10^^ cm. The magnetic field inside the hole corresponding to a 
half-flux quantum is Bo = {'kUc/ enB?) = 2 x 10^^ G, whereas the magnetic field in the junction is as high as 
Bj - BoL/d = 10-1 G. 

There are both nuclear spins and paramagnetic impurities in the spin bath. Consider first the nuclear spins; 
assuming all nuclei have spins, we find that in the bulk of the ring, within a penetration depth of the surface, there 
are Nr ~ 2ttRXlL x 10^'^ « 5 x 10^^ nuclear spins coupling to the ring current; and in the junction itself, a number 
more like Nj — (TrdZA/,) x 10^^ w 3 x 10^. Thus each ring nuclear spin couples to the SQUID with a diagonal coupling 

ujr ^ fJ-nBo ~ 2 X lO^i'^if; on the other hand for junction spins this coupling is Lu^j = fJ-nBj ^ 10^^ K. Notice we have 
ignored any coupling to substrate spins (assume, eg., the ring is in superHuid IIe-4!), which might have a much larger 
coupling to the current. 

At any temperature such that kT 3> ujI , uj^j , the typical polarisation of these spin baths will be \/N^, \J Nj re- 
spectively, giving a distribution of longitudinal bias energies with typical values E^ ^ aj| w 5 x lO'^'^K, and 
Eg ^ Loi^/ Nr ~ IQ^'^K acting on the tunneling flux coordinate $. If we now add paramagnetic impurities to the ring, 
with concentration ripm, and coupling w|m '^-^ 2 x IQ^LOr ~ 4 x lO^^*^ K to the current, this gives a typical longitudinal 

1/2 

bias energy E^"^ ~ npm x 0.2 K . This longitudinal term is obviously bigger than the nuclear contribution, unless the 
superconductor is very pure indeed! 

However, there is another much stronger transverse term, because each spin feels the dipolar fields from the other 
spins. This field is > \G (much higher near to the paramagnetic spins), and for the nuclear spins has an associated 
energy > lO^^X, which is ^ a;|!,a;j'. Physically, when the SQUID flips, the field on each nuclear spin hardly 
changes its direction, being dominated by the more slowly varying (but much stronger) nuclear dipolar field. For the 
paramagnetic spins the analogous coupling uj^^n > lO'^ times larger, which in the absence of nuclear fields would give 
an inter-paramagnetic "fiip-flop" rate V^™ ~ lO^ripm Hz, except that in pure samples these processes will themselves 
be blocked by the local dipolar coupling between the impurity and nearby nuclear spins (of strength LOp,-^ ^ IQ^'^K); 
this will happen once iipm ^ 10"^. 

We can write down an effective Hamiltonian for this system, valid over timescales considerably greater than A^^; 
we will use this later to analyse the effect of the spins on the SQUID dynamics (section 5.B). We will assume that 
A S> Vkk' (the only case of experimental interest); then we cauptxeat the internuclear dipolar fields as slowly- varying 
in time. The effective Hamiltonian can then be derived, to givellS: 

JV 

Heff(f^„) = {A„(<i>„)f+e-^S. "-^-^ + H.c.} + CffT, +Y.'^TJlal + c^^a^] (2.29) 

k=l 



11 



where wj! = uJr, or ojpm, depending on the spin, and has just been discussed; and microscopic analysis shows 

that \dk\ ^ ujl/flo, where flo is the Josephson plasma frequency (section 2.B). Notice that loI/uj-^ <C 1, which is 
the opposite limit considered to that for the giant spin! Notice further that these couplings are far less than 
(i?o, E^,EP"^ ^ Ac only because there are so many spins involved). Thus the spin bath is no longer "slaved" to the 
central system. In section 4 we see how this allows a mapping to an oscillator bath, coupled to f (ie., a spin-boson 
model). 

Finally we note that the external field also acts on the nuclear and paramagnetic spin dynamics, via the Zeeman 
coupling (which simply adds to oJb. This can help to suppress decoherence effects, by freezing out the spin bath 
dynamics- for more details see refJl3. 



H. General Canonical Models 



We now recall our assertion that almost any mesoscopic "central" system, coupled to its environment, may be 
described at low e nergie s by (2.2), with the environment being written as a sum of an oscillator bath term (2.3) and 
a spin bath term (( 2.12 ), or a higher spin generalisation). The simplest example is of a single central spin coupled to 
both oscillatpi' and spin baths. Such a model seems forbidding but in fact a fairly complete analysis has been given of 
its dynamicsB- we recall some of the results at the end of section 4. One can also consider a much more complicated 
model in which a macroscopic array of central spins {fj}, at positions {r^}, couples to both oscillator and spin baths. 
The effective Hamiltonian is then an obvious generalisation of what has gone before: 



jT^ COS 



N 



k 

N 



H.C. 



En 



N N 



Eu/Esfe^^"* •'^* + E 



k=l 



k=l 



E E N// + (s", V + h-c.)\ + ^ E ( ^ 
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J <i 



k=l k' = l 
1 



4- mqLo\x\ 



(2.30) 



where V(yi — rj)TfTj is a "high-energy" diagonal coupling between the various "central spin" systems. If we throw 
away the spin bath we get a set of 2-level systems coupling to aH|-escillator bath, of which the simplest example is the 
"PISCES" model (in which there are only two 2-level systemsESEJ) . 

Such models seem impossibly complicated, but actually one can solve for their dynamics in many important regimes! 
What is crucial is the separation of the 2 baths. Often (as with nuclear spins) their mutual interaction is very weak 
(and can be parametrised by a time Ti{T) which may be very long at low T); in this case this separation is a good 
one. If there are certain spin bath modes that interact strongly with the oscillators, then typically we can simply 
absorb these modes into an "augmented" oscillator bath by a canonical transformation. An obvious example arises 
with electronic spins in a metallic host (the Kondo or Kondo lattice problems); one rewrites the bath to include the 
"Kondo resonance" in the oscillator bath spectrum. 

A proof that one may do this in all cases seems rat her d ifficult- in any case the usefulness of these models tends 
to be established by ptheic-apphcation. Models like— ( 2.3C ) describe mesoscopic j-systems like coupled SQUIDs or 
coupled nanomagnetsL^LjGJ, Quantum Spin Glassea^j and low-T dipplar glassea22l, as well as coupled anisotropic 
coupled Kondo spins and Kondo lattices, coupled nuclear spin systems^, superconducting arrays, or coupled defects 
in solids. They are also useful for analysing purely theoretical questions about relaxation, dissipation, decoherence and 
quantum measurements in quantum systems- many oupstions remain unanswered . .h aving only been studied thoroughly 
in restricted models such as the spin-boson modeM'EJ or the PISCES modelHE^. We return to experimental and 
theoretical applications in section 5. 



III. AVERAGING OVER THE SPIN BATH 



To extract useful information from the low-energy canonical models, we must calculate their dynamical properties. 
Since we are typically not interested in the environment (one usually has little control over it), one performs a 
statistical average over the environment. This procedure is fraught with danger, because of "memory" effects in the 
environment, and because assumptions such as "self-averaging" in the environmental correlation functions may not 
strictly be valid. 
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In this section we show how the spin bath may be "integrated out" by means of 4 different statistical averages, 
each involving an integration over a particular variable. The end result is a description of the time evolution of the 
"reduced" density matrix for the central system- provided we can ignore memory effects-ip the environment. The 
starting point is no different from that involved in functional averaging over oscillatorglljil3; both begin with a path 
integral form for the propagator of the reduced central system density matrix, written as 

K{1,2) = dQ dQ' e-^/'*(^^''[01 - ^°[«'l).F[g,Q'] , (3.1) 

Jqi Jq[ 

where So[Q] is the free central system action, and J-[Q, Q'] is the famous "influence functional" El, defined in general 

by 

nQ,Q'] = l[{Uk{Q.t)uliQ',t)) , (3.2) 

k 

Here the unitary operator Uk{Q,t) describes the evolution of the k-th environmental mode, given that the central 
system follows the path Q(t) on its "outward" voyage, and Q'{t) on its "return" voyage; and !F[Q,Q'] acts as a 
weighting function, over different possible paths {Q{t),Q'{t')). 

For a central 2-level system, the paths Q{t), Q'{t) are simple (recall Fig. 2): 

2n 

Q(n)(s) = 1 - ^ [sgn{s - t2i-i) + sgn{t2i - s)] , (3.3) 

i=l 

where sgn(x) is the sign-function, and n is the number of transitions of the central system, occuring at times 
ti,t2, ■ ■ ■ ,t2n (for definiteness we assume trajectories starting and ending in the same state, and use the conven- 
tion that Q = ±1 corresponds to Tz = ±1). The goal is to find the central spin density matrix; in this article we give 
results for the "return probability" Pii{t) for the system to be in the same state | t) at time t as it was at i = 0. 



Using (3.3) this can be written as an "instanton expansion" over flips of the central spin (Appendix A): 

Pll(t) = ^(*A„)2("+™) / dt^... dt2n / dt[... rf4„-^[Q(n): (3-4) 

Further simplification arises if the environmental modes are uncoupled- then J-[Q, Q'] factorises, and we can write 

^[Q:Q'] = exp(— Q']) = exp(— i Q']), where the complex phase (j)k[Q,Q'] contains both real (reac- 

tive), and imaginary (damping) contributions. 

Now for an oscillator bath one simplifying feature is crucial, viz., the very weak coupling to each oscillator. This 
allows one to evaluate each up to 2nd order only in these couplings, in terms of a spectral function for 

the unperturbed oscillator dynamics (compare J{oj,T) ia-.(^)). Even though the paths Q{t) and Q'{t) may be 
complicated, the calculation of J-[Q, Q'] is often tractablellil. 

However because the coupling to each spin bath mode is not necessarily weak, it will in general strongly alter their 
dynamics, often slaving them to the motion of the central system. Thus we cannot start from the unperturbed spin bath 
dynamics- the problem is fundamentally non-perturbative in the {wfe}. However it is not intractable, because one can 
rather easily deal with the dominant longitudinal terms {wjl}. The other terms can then be dealt with perturbatively 
(and sometimes even non-perturbatively) . It is the separation of the effects of the various terms in the Hamiltonian 
which leads to not one, but 4 different averages. What is quite remarkable is that these averages can be evaluated 
analytically in most cases (section 4). 

We begin by explaining the 4 different averaging integrals required for a general spin bath. This is done pedagog- 
ically, by solving for 4 different limiting cases of the central spin Hamiltonian (sections III.A-III.D), each of which 
requires only one of the 4 averaging integrals. Then the general procedure (combining the 4 averages) is given in III.E. 
One reason for going through these averages one by one is that each corresponds to a different physical mechanism 
of decoherence- we return to this in_s£ction 5 (note that more detailed results for the 3 limiting cases discussed in 
sections III.B-III.D are given in refs.Lll3). 

A. Phase averaging: Topological decoherence 

Formally the case of pure topological decoherence applies to the following special case of H^ff, in which only 
non-diagonal terms are included: 
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TV 



CTT 



2Ao I f_ cos - « ^ otkUk ■ Sk] + H.c. 



(3.5) 



k=l 



Since lo'^ = lo 



^ = are zero, all the 2^ environmental states are degenerate, and there is no exchange of any energy 
between r and the {c^k}- The only thing that is exchanged is phase] the phase $0 of f becomes entangled with that of 
the {(Tfc}, during the transitions of r between |) and |), so that the initial and final states of the spin environment are 
different. Physically ( |3.5| ) would arise if the original high-energy Hamiltonian contained only "transverse" couplings 
(to {c^}), which only act while the central system is making transitions (in the case of a moving particle coupled to 
a spin bath, they would be " velocity couplings" , only acting when the particle is moving) . 

We wish to determine Pii{t) for this case. For pedagogical purposes let us begin by assuming that —iak is real, ie., 
Uk is pure imaginary; then we have added a pure environmental phase term to the free 2-level Hamiltonian. Then, 
writing —iak ~> c^k^ the formal solution to this problem can be written immediately asEl 



Pn(t)4 1 



(cos [4Aot cos ($0 



N 



k=l 



ctkrik ■ CTk 



(3.6) 



where the bracke|ts_(--) trace over the spin bath. By writing this as an instanton expansion over central spin flips (see 
App. A and refs.EH3), we transform it to a weighted integration over topological phase: 



00 „2 



dip 



m— — 00 



1 1 



1 f °° 

- 1+ E (~irp( 

m— — 00 



C08{2Ao{ip)t) 

™)e"™'J2„.(4A„ol , 



2n 12 2 



(3.7) 
(3.8) 



where Ao(ip) = 2AoC0S(p is a phase- dependent tunneling amplitude, J2m{z) is a Bessel function, and 



N 

F{m) = {We 

k=l 



N 



cos(2mQ;fc) . 



(3.9) 



k=l 



For small Uk we may approximate the product in (|3.9| ) as 



N 



(3.10) 



k=l 



Notice that A is just the mean number of environmental spins that are flipped, each time t flips. 

Clearly phase decoherence is important if A > 1, in which case F\{y) = Sufi + small corrections. Then, rather 
surprisingly, we get a universal form (shown in Fig. 6), in the intermediate coupling limit, for Pii{t): 



Piiit) 



1 + Jo(4A„i) 



/dip 
-^Pi^\t^^ — if) {intermediate coupling) 



(3.11) 



with a phase integration over the free central spin propagator P^^\t,^) [cf. eqtn. (AS)]. Thus random phases 



arise because successive flips of t cause, in general, a different topological phase to be accumulated by the spin 
environment. In fact, the universal behaviour comes from complete phase phase randomisatioroQ, so that all possible 
phases contribute equally to the answer! The final form shows decaying oscillations , with an envelope ~ t^^^^ at 
long times, which can also be understood by noting that the "zero phase" trajectories contributing to Pu constitute 
a fraction (2s)!/(2''s!)^ ~ s~^/^ of all possible trajectories, where s ~ Aot. 

In the strong coupling limit of this model, where the bath spins rotate adiabatically with the the central spin, 
one has ak 7r/2, so that F{m) = (-1)™ and Pii(i) = ^[1 + cos(4Aot cos $)], where $ = $0 + NTr/2, i.e., the 
Haldane/Kramers phase is now $, since the N bath spins are forced to rotate with f. 

The results for complex ak are given in Appendix A; the basic ideas behind them (and the techniques for their 
calculation) are simple elaborations of the above. 
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B. Average over longitudinal fields: Degeneracy blocking 



We now consider the effective Hamiltonian 



N 



(3.12) 



k=l 



To solve th is we will a ssume the model discussed previously, in which all cluster around a single central value lOq 
(cf. eqtns (2. 24), (2. 23), and Fig. 4). Since the bath now just acts as an extra longitudinal field, we are dealing with the 

trivial case of a biased two-level system, with bias energy {S^h + e), where e = '^^=1 ^\^k- "^^^ only question is how to 
average over the internal bias- this depends on whether we deal with a single central system, or a statistical ensemble 
of them (corresponding to either an average over many measurements on a single system, or a single measurement on 
a large number of non-interacting systems- interactions are discussed in section 5. A). 

For a single central spin, the dynamics of r in this model are completely trivial- one has -Pii(i) = [1 — 
(A^/ £^^) sin^ Et\, where i5 = -I- e, ie., resonant tunneling of an isolated spin in a longitudinal bias feld (com- 
pare (|A4|)). 

For an ensemble of central spins, we must average over the whole bias range. In what follows let us assume for 
definiteness a spin bath at some equilibrium temperature T = 1/(3] then the ensemble average is just a weighted 
average over bias: 



where Z{j3) is the appropriate partition function. 



-/3£ 



Pii(i) - / deW{e) 



1 - 



de W{e) 



-f3e 



{e + inY+AAl 



1 - cos(2tv/(e + e//)'+4A2) 



(3.13) 



(3.14) 



The physical interpretation is obviousB; only a very small fraction A{(^h) ^ '-^o/Eq of central spins in the ensemble in 
the "resonance window", ie.,, having total bias |(^// + e)| < ^o, can make transitions- this selects states with internal 



bias around e ~ —^h- All other states lack, 
tunneling- they are "degeneracy blocked" Q 



e near-degeneracy between initial and final energies required for resonant 
The resulting correlation function is then 



Pii(t) = 1 - 2Ai-^H) ■hk+ii'iAot) 



(3.15) 



A;=0 



where A{^) = AoW^(^) (except in the unphysical case where the polarisation group linew idths rj\/ < Aq. For the 
usual case where all polarisation groups overlap, and W{e) has the Gaussian form (2.23), one has A(^)/(27r)^/^ ~ 
(Aq/Eq) exp(— 2^^ /£'^) . It is not surprising to find that the spectral absorption function x"{'^) — -^"^ / dtPii{t), 
corresponding to (3.15), has the "BCS" form 



8Ao 



16A2 



(3.16) 



Finally, let us note that one may imagine a case where one has an ensemble of systems in which, altho ugh the bath 
state is not fixed, the polarisation group is known to be equal to M. In this case we must replace ( 3.13 ) by 



de GM(e) 



-I3e 



(3.17) 



in (3.14) for Pii(t), where Zm{P) is the partition function for the A/-th polarisation group. If we then recalculate 
x"(w) in the same way we find almost zero a bsorp tion unless \MuJo + £,h\ < Max{Ao,TM), where Tm is again the 
linewidth of the Af-th polarisation group (cf. (2.22). 
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C. Average over transverse fields: Orthogonality blocking 



Until now we h ave ignored the "transverse field" part "^fe" ' '^fe the diagonal term in the effective 



Hamiltonian (2.16). To study this let us consider again an effective Hamiltonian which has no non-diagonal terms 



apart from the "bare" tunneling, but having all diagonal terms: 

N N 

HeS = 2AoTx + TzUjI ^ /fc ' CTfc + ^ CJ^ TTlfc • CTfe , (3.18) 

fc=l fe=l 

We will assume ul S> i^;^, ie., that the transverse "orthogonality blocking" part of the diagonal interaction is much 
smaller than the longitudinal part. To make things as simple as possible we drop all degeneracy blocking effects, ie., 
we assume all w|| are equal (ie., = 0). It then follows that the spin bath spectrum is split by lOo into 2N "polarization 



groups" of degenerate lines, with C^'^^''^'^ degenerate states in polarisation group M (cf. eqtn. (2.21) 



At first glance it seems that the {dk} in (3.1J) simply act on the central spin r as an external field. However 



this is wrong; it ignores their role as dynamic quantum variables. The dynamics come because th e \(fk } can precess 



in the fields {"fk} acting on them, and these change each time t makes a transition (cf. eqtn. (2.19)). Quantum 



mechanically, the precession caused by is equivalent to saying that some bath spins are flipped when the central 
spin T flips (in general a different number of them during each transition of r) . To see this formally, recall that cw^ 

exists when the initial and final fields 'y^,^'' and j^^^ acting on dk are not exactly equal and opposite (cf. eqtn (2.19)) 
Defining the small a 
of dk are related by 



Defining the small angle Pk by cos2/3fe = ~lk^^ ■7fe^Vl7fel)ll7fe2)| (recall Fig. 3), we see that the initial and final states 



I 4) = Uk I dr) = e-^"^'^^ I din ■ (3.19) 

Suppose now the initial spin bath state belongs to polarisation group M . If, when r fiips, bath spins also fiip so 
that M — M, then since E^{M) = Eij_{—M), resonance is still preserved, a transition is possible- indeed r cannot 
flip at all unless there is a change in polarisation state of magnitude 2M. For this change in polarisation of 2M, at 
least M spins must flip; moreover, for resonant transitions to continue (incoherently), the bath polarisation state must 
change by ±2M each time f flips. 

Let us therefore deflne Pnit) as the correlation function Pii{t) restricted to systems for which the bath polarisation 
is M . For a thermal ensemble, 

N 

Pn{t;T)= ^ w{T,M)PM{t), (3.20) 

M=-N 

with a weighting w{T,M) = Z'^C^f'^^^^^^e"^" "/''^^ , where Z is the partition function. 



In Appendix A. 2 we calculate PM{t) (see Eq. ( A27 ) as a weighted average over an orthogonality variable x: 

(3.21) 



PM{t) = ^ cfcx e ^1 + cos[4:Aq.JM{2xy/K)t]^ 



POO 

= 2/ dcxe-'''p^"\t,AM{x)) , (3.22) 
Jo 

where Pl^\t, Am) is just the usual free 2-level correlator (Eq.( |Al[ )), but now with an an cc-dependent tunneling 
amplitude Aj\/(a;) = 2A^Jm{2x^/k); and k is the "orthogonality exponent", defined by 

N 

e'"^Y[cos(3k - 6^(1/2^1^^^', (3.23) 
fc=i 

The orthogonality blocking term Pk is analogous to the topological decoherence term . It is important to understand 
why we must introduce the average over x. Mathematically, it comes from the restriction to a single polarisation 
group (see Appendix A. 2). Physically, it corresponds to a phase average just like that in ( ^ ) (compare the Bessel 
functions), but now the phase is that accumulated between transitions of the central spin (rather than during these 
transitions, as in topological decoherence). This phase accumulates if cj^ ^ 0, because then the field •jk on the fc-th 
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bath spin does not exactly reverse when the central spin flips, and so this spin must start precessing in the new field. 
It is random simply because the waiting time between flips is a random variable, iur-the path integral. 

We shall not give full details for the dynamics of this limiting case (for which see refS), but just enough to understand 
the physics. First, note that the terms Puit) are easily verified to be incoherent, and so Pii{t) ~ fPM=o{t) + 
incoherent^ where / = ^J2/ttN . Even the small fraction / of systems in an ensemble having M — Q will only have 
coher ent dynamics if k <C 1. The easiest way to see this is to again calculate the spectral absorbtion function from 
(|3.2lD, to get: 



E 



2Aok1/2 ^ |Ji(2Ki/2a;^.)| 



(3.24) 



which leads, as k increases, to an ever-increasing number of square-root singularities in xX/=o(^)- -'^'-'^ only 
a single root xi - [(1 - w/2Ao)/«;]i/2 enters, and x'm^o{<^) ^ (7r//2AoK)e[(i-"/2^°'/'"l for lo < 2Ao, ie., a fairly 
sharp asymmetric peak at the resonant frequency of the free spin. For larger k the multiple peaks and tails mix, and 
Xm=o('^) shows no obvious peak- moreover, its spectral weight is shifted to ever-lower frequencies: 



PM=o{t) - 1/2[1 
Ae// = 2A$e-'' 



- cos(4Ae//t)] 



K < 1 



(3.25) 



P,{t) = l-4Al^.ft;' + 0{At^.ft*) 



1/4 



K > 1 (3.26) 



Notice that this reduction of the transition rate is much slower than the usual polaronic or Anderson/Hopfield 
orthogonality catastrophe, relevant to oscillator baths, which gives exponential suppression of Ae/ f for strong coupling. 
This is understood as follows. In oscillator bath models, band narrowing comes essentially without any bath transitions 
(most of the polaron "cloud" is in virtual high-frequency modes) - it is adiabatic. Here, however, roughly k spins flip 
each time S flips (the probability of r flips is K^e~'^ /r\, which peaks at r ~ k), even though we only consider PM=oit), 
i.e., even though AM = (just as many bath spins flip one way as the other). 

A further examination ftWhc correlation function shows that the structure of Pm=o {t) is exceedingly bizarre- it was 
described in detail in rcf.Lm. In section 4. A and Fig. 8 we will return to the physics of orthogonality blocking, but 
including other mechanisms as well (see also Fig. 7 below). 



D. Averaging over spin bath fluctuations 



The previous 3 averages assume no intrinsic spin bath dynamics- the bath acquired its dynamics from the central 
spin. Consider now a Hamiltonian 

N N N 

H,s = At, + Cr, + ^ 4 -^fc + E E Kk^'^k4 ; (3-27) 

k=l k=lk' = l 

in which we assume \Vkk'\ but arbitrary Vkk'/A. The addition of Vkk', to what would have been a simple 

degeneracy blocking Hamiltonian, gives the spin bath its dynamics, and causes 2 changes (as noted previously in our 
introductory presentation of the central spin model). First, a polarisation group M acquires an "intrinsic linewidth" 
To ^ VqN^^'^, where Vo is a typical value of \Vkk' \ for the N bath spins (of course normally Fq <C I'm, unless the 
{uil} happen to be extremely tightly bunched together). Second, the transverse part of V^f^ causes pairwise flipping 
amongst bath spins (eg., transitions | tfcife') ^ I ifetfe')); ^ I'^^e ^ NT.^^, where Tj"^ ~ Vq. This "spin diffusion" 
in the bath causes the internal bias e to fluctuate in time, inside the energy range of the polarisation group M, with 
a random walk correlation ([e(t) - e(t')]^) = A^|i - t'\, where A^ = f^T^^, for timescales K\t - t'\ < 1. 

It thus follows that for a given single central spin, with its surrounding spin bath in polarisation group M, the 
problem reduces to calculating the dynamics of a 2-level system in a longitudinal bias field (^ -|- e(t)), where ^ is 
the applied bias, and the internal field is e{t) = Mujo + 5e{t). The correlation properties of 5e{t) are those just 
described- our task is simply to functionally average over these fluctuations in the calculation of Pii(i). In doing 
this we will make the physically sensible assumption of "fast diffusion" of (5e(i), such that the time At it takes for 
the bias to diffuse across the "resonance window", of energy width A, satisfies At ^ 1/A. Then the system has no 
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time to tunnel coherently, but can only make an incoherent " Landau- Zener" transition. Since the bias changes by 
5e ~ 5lOo{N / {T2l^)Y^'^ in a time 1/A, this formally requires that 

< = f 2^2"^ {fast diffusion). (3.28) 

This problem is solved in Appendix A, by performing a weighted average over dynamic bias fluctuations, with the 
restriction that these only occur inside polarisation group M; the relevant average is 

X>t{t) e"^ -f '^^^ -f ''*2A:(ti-t2)c(ti)e(t2) . ^2 29) 

where 2K~^(ti — t2) = A^{\ti\ + |<2| — \ti — ^2!) is the correlator of the dynamic spin bath fluctuations (see Appendix 
A). One finds that -Pii(i) decays as a simple exponential -Pii(i) = e~*/^'^'^^\ where 



r-/ = 2vri/2^e-(«+*^'^°)'/r' (3.30) 



This result is easily understood- the bias fluctuations can cause the system to pass briefly through resonance (allowing a 
transition of the central spin), but only if the net static bias ^ + MlUo is not greater than the range F of the fluctuations. 
By comparing with the case of pure degeneracy blocking we see that the important role of the bath dynamics is (i) 
to unblock the central spin dynamics, by helping it to find resonance, now over an energy window of width F (instead 
of Ao) around zero bias (recall Fig. 4), and (ii) to change the central spin dynamics from coherent to incoherent 

tunneling. 

Note that in a model like ( |3.27 ) we have eliminated bath fluctuations between different polarisation groups. The 



basic assumption is that any "Ti processes" in the intrinsic bath dynamics, which could change M in the absence of the 
central spin, are very slow (At low T, Ti for nuclear or paramagnetic spins does become extremely long). However this 
is not always realistic- we return to this point in section 5. A. If Ti is short, one must make a dynamical average over 
2 kinds of fluctuation, usually with quite different time correlation, viz., the intra-polarisation group fluctuations, 
described by ( |3. . 2(il ) .and with correlation time NT2^, and the inter-polarisation group fluctuations, occuring on a 
timescale Ti (cfOlJ). 

E. Averaging over the Spin Bath: General Results 

We now turn to th e problem of averaging over the spin bath for the general form of the Central Spin Hamiltonian 



given in eqtn. (2.16). This can be given in the form of a marvellously simple prescription- one simply applies the 
4 averages we have just seen, to the problem of a simple biased 2-level system! We begin by giving the explicit 
prescription (whose proof is given in Appendix B), and make a few comments on it. 

The prescription begins with the following 4 averages (all of which we have seen in the preceding 4 sub-sections): 

(a) A "topological phase average" V i^A'(i^) / ?e*2''(*-¥') ; (3.31) 

_ J 27r 



z^— — 00 



r 2 

(b) An " orthogonality average" 2 / dxxe~^ ; (3.32) 

Jo 

(c) A "bias average" j de GA/(e)^-^ OR J de W{e)^^Y. ' (^.33) 



(d) A "bath fluctuation average" j Ve{t) / '^'i / . (3 34) 

As before, we assume a thermal distribution over spin bath biases, with a corresponding partition function and Z{0). 
All averages are normalized to unity. The weighting function Fyiy) = e"^"^ " in (3.31) is a generalisation of that in 
( |3.10 ), to allow for arbitrary directions of the unit vector n^; we now define 
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^ = 2 E ■""•■'(l - K)') ' ^' = 2 E "'K)' ' (3-35) 

fe k 

Now, suppose we want to calculate Pii(i). The prescription is fairly obvious in the light of the results given above 
for the 4 limiting cases. One follows the following steps: 

(i) Begin with the quantity 

; Am x) ; e) - 1 - ^|4^ sin^ {Em (<^, x)t) , (3.36) 

which is just the free central spin correlator (cf. ( |A4[ )) having tunneling matrix element AmC^^, x)^ and in an "internal 
field" bias e. The energy splitting Em is given by E'^,j(ip, x) = A|^(9?, x) + e^, and the matrix element Aj\/ is 

^m[v.x) = 2Ao|cos(^)JA/(2a;/7)| , (3.37) 



A if A ^ K (topological decoherence regime) 
K if K ^ A (orthogonality blocking regime) 



(3.38) 



We defined k previously (eqtn. (3.23). We will not give results for the case k, '-^ \\ they are extremely complex, do 
not appear to add new physics, and seem unlikely to be realised in practice. 

(ii) Now carry out the averages over topological phase [Eq. ( 3.31 )] and orthogonality [Eq. ( |3.32| )], to give an 
expression PM{t, e) describing the central spin dynamics in a bias e, coming from a bath in polarisation state M: 



PM{t;e) 



2 r dxxe--" V F^,{v) f 



(*-<<=) 



1/ — — 00 



Ai/(y>,x) 



S\T?{EM{f,x)t) 



(3.39) 



where the weighting function is Fxi{v) ~ e over winding number v (recall eqtn. ( 3.10D ); 



(iii) Then, carry out the bias average [Eq. (p^)]. We wiU assume in the following for definiteness an ensemble 
average over all polarisation groups, thereby ensuring a summation over M, to give 



-I3e 



N 



Pn{t;T) = l- I deW{e)^^ ^ {1 - PM{t,e - McOo)) ; 



(3.40) 



M=~N 



This result summarizes the central spin dynamics in the case where the spin bath has no dynamics of its own, and 
only acquires dynamics through its interaction with the central system. In some cases there will be no intrinsic bath 
dynamics, and this will be the final answer. If we wish to apply the theory to a single central system, o r for some 
reason we can fix the polarisation group to be a definite value M, then we drop the summation over M in ( [3.4C| ), and 
replace W{e) by GM{e). 



(iv) When the interaction ter™, Vkk' plays a role, we apply the 4th average (3.34) to (3.40), as described in Appendix 
A (cf. also the discussion in ref.cl). This gives the completely incoherent form 



/oo 
xdxe--' J2 / ;^Fv(;.)e^2''(*-^)[l 



~t/TM{x,ip)'\ 



(3.41) 



where the relaxation rate Tj^^{x, ip) is given by 



dse^' 



(3.42) 



with Am(2^, </?) given by (3.37), and where G^(e) is a Gaussian of width V = /ic^o- Since F 3> ^ we have also F 3> A, 
and so we get 
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Tm W = > (3-43) 

This result is the most general one for the dynamics of the central spin, if all 4 bath averages are included- it is 
generally valid, with only the single restriction that the diffusion of the fluctuating bath bias in energy space be fast 
(cf. eqtn. ( |3.28 ). In the absence of such fluctuations we go back to (3.40). 



In essentially all physically realistic cases the different polarisation groups strongly overlap. It is then simplerE 



to transform the sum over M in (3.4C) or (B.42) into an integral over energy bias ^ , using the change of variableja 
X^Af ~^ J d(,/2ujo, and then integrate over ^. One way to do this (using steepest descents) was detailed in ref.H 
(compare eqtns (4.41)-(4.47) in that paper). For exact answers one can use the identity 

/ xdxe-''\jlj{2x^)cos^ (t) = lM{'2j)e-^'^ cos^ cf) (3.44) 
Jo 



to evaluate either ( 3.40| ) or (3.41) 



In the next section we will evaluate Pii(t) and its Fourier transform for a number of different parameter regimes. 
But even before doing the integrals, the qualitative behaviour is obvious. Relaxation is only occurring for central 
spins which happen to be within a bias of exact resonance. The width of this "resonance window" is coming 
from the energy which the bath spins can provide to the central spin, by flipp ing up to ^ 7 bath spins; hence ~ 7Wo 



(formally this is obvious from the properties of the Bessel functions in ( 3.37 ) and ( 3.44 ), which fall off very fast once 
M > 7). We show graphically the relaxation of different groups in Fig. 7(a); again one sees how only groups with 
M < 7 relax quickly. The T2 bath fluctuations help this process by bringing the a central spin in polarisation group 
M to its resonance window (of width Am)- Only transitions of systems having M = can show (partial) coherence; 
all transitions with M are essentially incoherent. Note that the resonance window will not be visible in a resonant 
absorption experiment (Fig. 7(b)); higher M groups contribute only a very low frequency contribution to this. This 
nicely demontrates that one is very far from any linear-response regime in the present system (so that, e.g., the 
fluctuation-dissipation theorem is somewhat irrelevant here). 

Without the bath, transitions of the central spin would be coherent, but over the far smaller resonance window of 
width ^ Aq. If we only had T2 bath fluctuations, but 7 <C 1 (ic., no bath spins flipped during the transitions of the 
central spin), then we would again get incoherent relaxation, but this time with ~ rA/=o (the width of the M = 
polarisation group). One can also imagine a situation in which Ti is very short, so that all polarisation groups are 
involved in the relaxation, and the resonance window is just the distribution Vt^(C), with — Eo — N^^^coo- 

These results thus tell us that in the presence of a spin bath, any ensemble of cgptral spins, initially spread over a 
range of biases, will start relaxing by digging a "hole" of width around zero biasBH. This hole reflects the intrinsic 
central spin dynamics (ie., it is not being produced by interaction with some external resonant signal- it should not be 
confused in any way with the "spectral hole-burning"^^d£me by experime nters working on glasses or in optics, using an 



external source). As discussed in a number of paperslo, evaluation of (3.41), using either steepest descents or other 



means, shows that the system under most conditions relaxes incoherently with a relaxation rate (after summing over 
all polarisation groups, and doing the orthogonality and phase integrals) given approximately as a function of bias ^ 

by 

TT^'-^l 

All of this is in complete contrast to how inelastic tunneling works in the presence of an oscillator bath0; there 
the relaxation rate typically increases as one moves away from resonance, usually as a power in bias (t~^(^) ~ ^ 
for diagonal coupling to phonons, ~ for non-diagonal coupling to phonons, and ^20-1 diagonal coupling 
to Ohmic baths like electrons via a dimensionless coupling a). Thus one does not expect hole-digging for oscillator 
bath-mediated quantum relaxation, except over a very narrow region of width ^ Aq. 

Finally, we note that one may also give a formal prescription for the case where some central spin couples simulta- 
neously to an oscillator bath and a spin bath. We do not give the details here- they are discussed fairly exhaustively 
(along with the results for the central spin dynamics) in ref-EI (see also end of section 4). 



IV. DYNAMICS OF THE CENTRAL SPIN 



Given the large number of parameters entering into the 4 averages just described, we seeJittle point in an exhaustive 
description of Pii{t) over the whole parameter domain (for more extensive results see refs.Q'u). Instead we concentrate 
on 3 points. First, we show how in the strong coupling regime, coherence is destroyed, leaving incoherent quantum 
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relaxation; this regime applies to almost all mesoscopic or nanoscopic magnetic systems, because of their coupling to 
nuclear spins and to paramagnetic impurities. Second, we discuss the physics of the weak coupling regime (applicable 
to, eg., SQUIDs), and how in one limit of this regime one may formally map the spin bath to an oscillator bath. 
Finally, and very briefly, we comment on the results obtained when one couples simultaneously to a spin bath and an 
oscillator bath. 



A. Strong coupling regime 

As already explained, the strong-coupling regime is defined by the condition wji and/or > A,,. This condition 
applies to virtually all situations in which the couplings are hyperfine ones to nuclear spins, or exchange couplings to 
paramagnetic spins; and also when one has dipolar couplings to paramagnetic impurities or defects. 

Almost all interesting physical examples in this regime fall either into the catagory of "strong orthogonality blocking" 
(when K ^ X') or strong "phase decoherence" (when X ^ k). In both cases the central system makes transitions 
accompanied by flips in the bath spins- so that even if the isolated central system is not in resonance, it can "find 
resonance" by using the flipped bath spins to make up the energy difference. If the couplings are such that roughly 
K bath spins flip, the range of energy bias over which transtions can occur is extended to roughly 2klUo- The central 
system is helped in this task by the fluctuations in bath bias caused by the interspin interactions Vkk' ■ 

In what follows we concentrate on the physics of decoherence in this regime, with an eye to the physics of "qubits" 
and of " macroscopic quantum coherence" . We also look at the form of the relaxation. We begin by explaining the 
results without the bath fluctuations, and then show what happens on adding these. 



(i) Results without bath fluctuations: In this case we must evaluate ( 3.40 ) and ( 3.39 ), suppressing either the 
orthogonality average or the topological average. In what follows we look at each case in turn, focussing particularly 
on the M — polarisation group contribution. 

(a) Orthogonality Blocked regime (k ^ X). In this regime only the orthogonality average 2 J xdxe~^ , and the 
average over bias e, are relevant- the phase average is approximated by a delta-function. The presence of the x- 
dependent transition matrix element Am (a;) — 2Ao J7\/(2x^/k) means that polarisation groups with M ^ k play a 
dominant role. Suppose however that we are interested in any coherent dynamics of the central spin- what will be 
found? It is obvious that transitions with finite M will be essentially incoherent, so we concentrate on central spins 
for which M = 0. Thus we simply integrate PM{t, e) over e, in the weighted bias average, to get 

PiY="(t) ^l-4A f dxxe^'^' I Jo{2xV^) \ ^ Jsfe+i [^K \ .h{2x^) \ t] (4.1) 

^ 1 n 



fe=0 



1 - 2 / dxxe-'^^ 2A{x) V Jafc+i [2Ao{x)t] , (4.2) 



fc=0 

where the x-dependent spectral weight is A{x) = A\Jo{2x^/K)\. Notice we have just done an "orthogonality average" 
over a "biased averaged" expression for the free system with x-dependent tunneling frequency Ao(a;). A Fourier 
transform to frequency space (which is essentially a picture of the relaxation rate, as a function of energy bias ^ for 
this system) gives the absorption spectrum 

XM=oi^) = lj dxxe-^\A{x) ^^^ i^AgtilV^ ^^^^ " ^ ' ^"^"^^ ' ^^'^^ 

Fig. 8 shows some representative plots for this "coherent" part of x"(u;); it is in fact almost completely incoherent, 
with total spectral weight 



(4.4) 



{dijj/2TT)^'{uj) =2A dxxe-"" I Jo{2y/lix) \ 
J J 

_ 2r(3/4) A 

- -3/2 ; 



a result which is very accurate even for k ~ 0.02. Note that the shape of x"(w) will change once we include all 
other (incoherently relaxing) polarization sectors M 7^ 0, and its total weight increases - in fact the total weight is 
^ Ak^>^ for large k, since ~ k^/^ different polarization sectors contribute. The absorption in x"(^) from these higher 
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M groups will be concentrated at frequencies ^ Aq (compare also Fig 7(b)). Note however that relaxation itself (not 
described by the linear response function x"(a;)) will be spread incoherently over a frequency range 2k,uJo, ie., 
the "hole-digging" in the relaxation occurs over a window of width IklOq. 

(b) Phase decoherence regime (X^ k). Let us now suppose the transverse field terms {i-i^^} are negligible compared 
to the {a/c}. We then deal with an effective Hamiltonian 



{N ~j AT 

f_ cos - « X! ■ ^''l + \ +'^^^'^l^k ■ (4-5) 
fc=l J fc=l 

Since ojj! ^ Ao by assumption, energy conservation requires that environmental spins flip with the central spin, just 
as in our discussion of pure orthogonality blocking. Thus in this case we must also keep th e ort hogonality average, 
to enforce this constraint, ie., we must perform the full average embodied in eqtns. (3.39),(3.4C). The full answer, 
including both the real and imaginary parts of ak, is rather complicated, and is presented in Appendix B. Here we 
will consider the more transparent answer one gets when is purely imaginary and adds to directly as a random 
variable to the central spin phase. 

Let us again start with only the M ~ contribution to Pii(i). Then we have 

Pit^it) = / deWie)^Poit,e) ; (4.6) 

with Po(^i given by ( [3.40 ) with M — 0. We may now carry out the integration in (4.6), assuming that W{e) is given 
by the usual Gaussian form ( 2.23 ), to get 



<=°W = l-2 / dxxe-^' ^^'(^) / |^e^'''^*"'^^2A(v.,x)5]j2fc+i[2Ao(^,x)t] 

with A{ip,x) — Acosip Jo(2a;\/A). The corresponding absorption x"(ijj) is 



(4.7) 



l^ — —00 



(4.8) 



It is possible to write analytic expressions starting from (4^), but in this somewhat pedagogical presentation we 
simply discuss the case when ^ = 0, i.e., zero degeneracy blocking, when LUk = ujo for all nuclei. The integration over 

bias is then absent (since W{e) is now just a set of (5-function peaks, ie., W{e) J2m=-n ^'i^^^^^^^^i^ ~ ^'^^o)), 
and we get 

Pn (i) = 5] u;(r, M)Pm {t) ; w{T, M) = / Z[I3) , (4.9) 

M 

where Pm(^) now describes the dynamics in zero bias; it is given by exactly the same weighted average over phase as 
in (U): 

PAiit) = f dxxe-''' ^ FA'(m) /"^e'2'"(*-'^) {l + cos[4AoiJM(2a;\/A)cos</?]} (4.10) 

m—-—OQ 

= Jdxxe-"ll+ f2 {-irFy{m)e'^"'''J2m[4AotJM{2xVX)] \ , (4.11) 

This can be interpreted either as an orthogonality-blocked expression, with frequency scale AMi^,x) = 
2AoCos{tp)JMi2x^/X) which is then averaged over ip, to give phase randomisation; or as an integration J dx over 
an already topologically decohered function having frequency scale Am{x) — 2AoJm i^xy/X) . It is intuitively obvious 
(and easily demonstrated) that only Po{t) may behave coherently, with a fractional weight ~ y^2/TTN in an ensemble. 

There are various interesting cases of ( |4.11 ) for M — 0. If A = (i.e., fik is parallel to z), then we go back to 
pure topological decoherence - the projection operator then commutes with the cosine operator. On the other hand 
if A' = 0, we have pure orthogonality blocking as stated earlier, and in fact when A' = 0, the parameter A plays the 
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role of K in (3.21). Notice that whereas the case A 
it pertains to the model in Eqs. 



can only occur accidentally, A' = is quite common - indeed 



(12.261) . 

We really begin to see the analogy between orthogonality blocking and topological decoherence when A, A' 1; 
just as with pure topological decoherence, Fy (m) collapses to a Kronecker delta, and we get the universal projected 
topological decoherence form: 



dxxe 



1 + Jo[4Ao^Jo(2a;^/A)] 



(4.12) 



2 



dxxe 



[16A2j2(2a;VA)-cj2]i/2 



7?(4Ao I Jo(2xVA) I -u) 



(4.13) 



which generalizes the result of ( 3.11 ) for pure topological decoherence. Eg. ( 4. 12 ) should be compared to ( 3.21 ). 

We show in Fig. 9 some results for x"('^) for selected values of A. The results are startling; even a very small value 
of A significantly washes out pure topological decoherence; but for any large value of A', we never get back the pure 
orthogonality blocking spectrum. 



The results in the case where is real, and n/j is along the af-direction (see Eq. (2.28)) are obtained by simply 
converting the Bessel functions Jm to Bessel functions of imaginary argument Im- 

Finally, note again that the above discussion of the AI ~ Q polarisation group is irrelevant to the real experimental 
lineshape- an evaluaticm of the full expression (3.40), summing over all M, for the A, A' ^ 1 regime, simply gives 
incoherent rclaxationclEj, spread over a frequency range ~ Acjq. 



(ii) Including Bath Fluctuations: The modifi-catioa of the above results, occasioned by the intrinsic spin bath 
fluctuations, was given in detail in Prokof'ev and Stampa. The fluctuations in bias allow the central system to cycle 
rapidly through the whole range of biases within a given polarisation group (transitions between different polarisation 
groups can occur through Ti processes- usually much slower). Here we simply recall the main result, which is ob tained 
by summing the relaxation forms from each polarisation group in an ensemble (cf. eqtns. (3.40) and (3.41)), and 
assuming that the spin bath Ti is longer than all experimental times scales. 

For a single central system, coupled to a spin bath in equilibrium at temperature kT ^ tJo, one finds that after 
an initial short-time transient, the relaxation is roughly logarithmic over a very long period; in fact one finds for the 
strong coupling regime that 



Hi/To) 



2ttN 



In 



{t < ic) 



(4.14) 



for times t ^ tc, where t 
group (copppare eqtn 



To{2N/e'^jy^^ , and Tq = 2A2/7r^/2p jg ^j^g relaxation time of the AI — polarisation 
(3.43)); thus tc is extremely long! For t ^tc the system settles down to a rather different 



behavioum. This logarithmic behaviour can be roughly understooda as coming from a. distribution of barrier heights, 
for the different polarisation groups, which are then summed over- as discussed in ref.u, section 4.3(b), the final result 
looks basically the same as that shown in Fig. 7(a). The fastest relaxation comes from the those polarisation groups 
in the "resonance window" (recall the discussion at t he en d of section 3). 

Two cautionary notes are in order here. First, (4.14) applies to a single relaxing system- but in the case of 
nanomagnetic systems, all experiments until now have been done on large numbers of nanomagnets, coupled together 
via long-range dipolar forces, which drastically changes the relaxation (see section 5. A below). Second, ( 4.14 ) should 
not be applied uncritically to experiments, even on single quantum systems. This is because in a real experiment 
there will also be (i) couplings to oscillator baths, and (ii) the relaxation -jiaill change once t >Ti. In superconductors 
or metals, electronic oscillator baths will often dominate the relaxatior£3, even at short times (for their effect on 
coherence, see section 5.C below). Even in insulating systems, one eventually expects the coupling to phonons to take 
over at long timesu, since this causes expo nent ial decay- even if very slow, this will eventually become faster than the 
spin bath-mediated logarithmic decay in (4.14). 



B. The weak couphng regime; relation to the oscillator bath 

A question of considerable theoretical (and practical) interest is the transition to the weak coupling regime, where 
the perturbation on the central system dynamics by a single bath spin is small (even though the net effect of all bath 
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spins, measured by parameters like A or k, may still be large if N is very large). The weak-coupling regime is thus 
defined by the condition ujk ^ Aq, ie., both and wj! are <C Aq. Again, a variety of cases is possible depending on 

how large are ratios like Vkk' /Aq and w^/wjl, or parameters like A and k. In the following we will not be exhaustive, 
but simply consider two theoretically interesting cases, in which Vkk' is assumed negligible, and we look at the limiting 
behaviour arising when cither w^/wjl <C 1 or cj^/cjj! ^ 1. We will also assume TV 1, otherwise the problem is 
trivial (the bath has little effect at all). 

There are 2 ways to solve for the dynamics in this regime. One is to use the averaging already developed above- 
this simplifies considerably in the weak- coupling regime. The other is to map the problem onto an oscillator bath 
one, and then use standard techniques to solve this. We demonstrate the 2 methods by solving one problem with 
each. 

(i) Longitudinally dominated ca se ( uj^/luI <^ 1): We assume the same Hamiltonian as in the discussion of the 

phase decoherence regime (eqtn. (Ej)), but now we can drop the orthogonality average- the projection to a single 

— II 

polarisation group is not required since wj! <C Aq. Again, for simplicity we consider the case when ak is imaginary. 



We shall solve this using the techniques previously developed; we let a; = in ( |3.37D , and hence use a matrix element 
A(ip) = 2AoCOS(/3 (which is independent of M). All polarisation groups overlap, and so we simply average over bias 
and topological phase: 

OO /"J ^ 

= 1- J2 F^{m) £e^'^^^-'^hA{^)J2j2k+i[2Ao{v)t], (4.15) 

m=-oo k=a 

with Ao((p) = 2AoC0S(/5 as before, and A{ip) — Acosip. This gives an absorption form 

X'\u;) = - £ F,{m) f ^e^^^C^-^) "^f^ ^ .,/, ^(W4A„- | cos^ |) , (4.16) 

^ J 27r [(w/4Ao)^ -cos^]!/^ 



which for large A simplifies to 

„, , 2A f dip cosmos / A I IN 

= — / lAfK\2 2 11/2 '? COS(p , 4.17 

uj J 27r [(w/4Ao) — cos^ (^J^/^ 

and can be expressed in terms of Elliptic functions. 

Notice that since uok <C Ao, the number N of environmental spins must be very large to have a noticeable effect, i.e., 
for A to be appreciable. Thus if A ^ ^OLk ^ since ak ~ ujk/^o, we have N ^ (ilo/ujk)'^ ^ (i^o/Ao)^. This not only 
implies that /i — N^^'^SuJk/^o 3> 1 (ie., very strongly overlapping polarisation groups) but also a Gaussian half- width 
N^/'^ujo ^ Ao, so that the internal bias e ^ Aq. The reason why the two mechanisms (topological decoherence and 
degeneracy blocking) are so easily combined is just because this bias is produced by all of the environmental spins, 
whereas only a few of them are actually flipped and their contribution to the bias field is small. 

(ii) Transverse dominated case (w^/tjl 3> 1): Consider now the case described by 

N 

Hcs = Aof, + Y^^ujidl + t.ujI&I) , (4.18) 

k=l 

with uj\lbjj: l.|--One can map this problem to the spin-boson problem, in a way similar to that given in the paper 
by Caldeira et al.E£l. We begin by noting that the time evolution operator for the fc-th bath spin, under the influence 
of a central spin moving along a path (5(„)(i), is just 

J/fe(Q(n),t)=T,exp|-j^ ds[u;ia'k + Q(n)is)Lula^k]'j , (4.19) 

Weak coupling means we can expand the time-ordered exponent to second order in wjl; completing the average 
in Eq. (3.4), and exponentiating the answer, one derives the influence functional for this problemlia. The "high 



temperature" result kT ^ cj^ is readily found 
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[Q(n)(s) 



is)] 



(4.20) 



The evaluation of this depends on the distribution of coupUngs. In the "strong decoherence" case where the bath 
states are spread over an energy range Eo (defined as usual by E'^ — ^J^ki^l)'^) which is large (ie., Eo 3> lu^, Aq), 
then J-^[Q, Q'] is simply approximated by its form for a single kink - anti-kink trajectory. With kink at i = and 
anti-kink at t2, one only has a contribution if uj'^{t2 — ti) ^ 1, which gives T[Q, Q'] = e^3^°(*^~*i^ , ie., a decay on 
a time scale E~^. It then follows that kink/ anti-kink transitions are bound in closed pairs, with phase correlations 



decaying over a time <C A^, ^, and the leading terms sum to give exponential relaxation, ie., Pii{t) = 1/2(1 
with a relaxation rate 



2Ao 



dte 



(4.21) 



On the other hand if Eo ^ Aq, we will get coherent oscillations of the system over a time scale ^ Ao/E^, ie., 
over roughly A^/Eg oscillation periods. We will use this result below, in discussing the possibility of seeing true 
mesoscopic quantum coherence effects in SQUIDs and nanomagnets (section 5.C). 

It is of course obvious from the original discussion of Feynman and Vernon that in the weak-coupling regime, whece 
an expansion to 2nd order in all bath couplings is sufficient, a mapping to an oscillator-hath .model must be possibleEj. 
The method given above is of course not the only one- others are explained in rcfs.Eilcl'HE3. 

Central Spin coupled to oscillator and spin baths: For reference we also briefiy note results for the case where 
a central spin is coupled simultaneously to a spin bath and an oscillator bath; the Hamiltonian is that in eqtn. (2.30), 
but without the sum over j, (and of course dropping the coupling term V{ri — rj) between different central spins). 
Some of the dynamic properties of this model were studied in ref.B- a complete study is still lacking. Quite generally 
one expects that at very short times the dynamics will be controlled by the spin bath, but at longer times incoherent 
oscillator bath- mediated transitions will take over. However these transitions still contain spin bath effects, via an 
integration over spin bath bias distribution, to give a result 

Pn{t,0=AZ-\p) f deM^(e-Oe-^(^-«){/me) + [l-/(r,e)]e-*/-(^'^)} 



(4.22) 

where we assume that the ensemble of central spina is in equilibrium with a bath of oscillators at temperature T = 1//3, 
in external bias ^ (for more general cases see ref.EI). In this case /(e) = e~'^'/2 cosh(/3e), and r(e;T) is the oscillator 
bath- mediated relaxation time; A is a " renormalisation" constant 0{1) coming from the spin bath dynamics. Most 
features of the results are obvious from ( 4.22 ). Thus, the oscillator bath unblocks transitions for a central spin way off 
resonance (ie., having ^ ^ Eg); these spins relax much as they would without the spin bath. When £^/leqEo there is a 
wide range of relaxation times, and the physical relaxation depends in a complex way on the oscillator bath spectrum 
as well as the bias. In any experiment one would see a crossover between purely spin-bath mediated relaxation at 
short times and this more complex behaviour- this crossover has not yet been studied in any detail. 



V. PHYSICAL APPLICATION OF SPIN BATH MODELS 



In this section we briefly review some recent work, mainly experimental, in which interactions with a spin bath play 
a role. We begin with magnetic systems where such interactions are strong, and where clear evidence exists for their 
influence on tunneling phenomena. We then discuss superconducting and normal systems. Finally, we discuss the 
difficult and controversial problem of decoherence and the mechanisms which govern it in nature. In this discussion 
the results of the Central Spin model are crucial- they describe the effect of a spin bath on a "qubit", or on a SQUID, 
or a magnetic macromolecule, which is trying to show coherent oscillations. 



A. Magnetic Systems 

(i) Magnetic Solitons: Magnetically ordered systems support a wide variety of soliton excitations, depending on 
the symmetry of the order parameter. These couple to various environmental excitations, which strongly affect itligir 
dynamicj&q These inoiwde both lineaB-and non-linear couplings to "quasiparticle" excitations such as magnonsriFi, 
phonon^Sl, electrons^, and photongia, all oscillator baths. More serious effects come from any localized modes 
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coupling to the soiiton, most notably paramagnetic impurities and nuclear spins. Both of these can be understood 
theoretically using a model in which a moving particle couples to a spin bathE3, as in section 2.C. 

Most experimental work in this area phaaJooked at domain wall tunneling in ferromagnets. Early experiments looked 
at the dynamics of multi-wall systema£3'Ej. More recentljj attention has focussed on the tunneling of single domain 
walls, whose position can be monitored in various waya^JED. These walls propagate along magnetic nanowires (of 
diameter 300-600 A)and are thus mesoscopic objects (the material of choice is usually Ni, for which the domain wall 
thickness Ap-|~ 700A). Although some predictions of the tunneling. iJaeory have been verified (eg., the square root 
dependencdij of the tunneling exponent onj-the pinning field He', cfou), the crossover temperatures to tunneling are 
usually much higher than calculated. valuesE^. A clue to the reason for this may be found in the microwave resonance 
experiments of Hong & GiordancO, which show extremely broad resonances, even down to 1.5 K. The coupling 
to phonons or electrons is far too small to give this; however. Oxygen impurities, which pill act as paramagnetic 
impurities, would act as a strong time-varying potential on any wall, and broaden the lineE3. It is clear that much 
more experimental work will be necessary to really understand the quantum dynamics of magnetic domain walls. 

(ii) Magnetic Macromolecules: On the other hand, a number of spectacular quantum effects have unquestionably 
been seen in crystalline arrays of magnetic macromolecules, in one of the most important developments in magnetism 
in recent years. Experiments on tunneling phenomena have had striking success in the "Mn-12" and "Fe-8" molecules 
(each of which behaves at low T as a "giant spin" of spin 10). There have also been highly-publicised experiments on 
the large "ferritin" molecule, which claim the observation of "macroscopic quantum coherence"; these are discussed 
in section 5.C below. 

The early molecular work demonstrated resonept tunneling (coming when snia_LeV[els of each molecule are brought 
into resonance), at relatively high temperaturea£3. More recent experimenta^J'L^CJ have gone into the quantum 
regime, in which only the 2 lowest levels of each giant spin are occupied. The classical-quantum crossover is clearest 
in the case of the Fe-8 system (which is very conveniently described by the easy axis/easy plane model discussed in 
section 2.F); below a temperature Tc ~ 0.4 K, the dynamics is completely independent of T (the energy gap to the 
next spin level in |tiiis, system is ~ 5K). In Mn-12 there are effects arising from "rogue" molecules which make things 
more complicatedHEa. At such low temperatures phonons are utterly irrelevant (experiments have now been pursued 
to r 30 mK), and the only dynamic environment left is the nuclear spin bath. The system thus seems at first to be 
an ideal realisation of the central spin model discussed in sections 3 and 4. The main complication is that experiments 
are done on a crystalline array of molecules, which interact via strong magnetic dipolar interactions. However, once 
one knows that the dynamics of a single central spin is going to be incoherent relaxation (cf. sections 3.E, 4. A), these 
interactions arc rather easy to deal withEj. The system is then described by an obvious generalisation of the central 
spin Hamiltonian, viz. 



in zero applied field (compare eqtn. (2.30), without the oscillators). Here H^^ is the central spin Hamiltonian ( 2.16| ) 
for the molecule at lattice position r^, and V{Ti — Vj) is the magnetic dipolar coupling between molecules i and j. 

To solve for the dynamics of such an interacting array of systems, one begina^j by defining a distribution function 
Pq,(^, r, t) for a molecule at position r to be in a longitudinal bias ^, with polarisation = a = ±1. It is then trivial 
to write down a BBGKY-like hierarchy of kinetic equations for Pa{£,,Y,t) and its multimolccular generalisations 
P(2)(l,2) = Pi?!a2(6,6;n,r2;t), and P^^\l,2,?,), etc., of which the first member is 

Pc (e, ^1 = -rN\0 [Pa ^1 - P-c T^] 

r. r') - p5 (e - aa'V{r- r '), f ; r, r') 



E 



df' f d£,' 



^0 J mie) 



(5.2) 



in which relaxation is driven by the nuclear spin-mediated relaxation rate Tj^^{S,) ~ (A^/f )e 1^1/^° (cf eqtn ( 3.45 ), 
in conjunction with the dipolar interactions. Under general conditions we must also solve for P^^-* in terms of 
P'^', etc.; but if the initial experimental state is either polarised or annealed then-.jES'^^ factorises at t = (ie., 
P(2)(l, 2) = P(1)P(2)), and can be solved. This led to the following predictionsBEl: 

i) Relaxation should only occur for molecules having |^| < ^o; consequently a "hole" rapidly appears in M(^,t) = 
P+ (C 5 ^) ^ P- (S,, t) with time, with initial width determined entirely by the nuclear dynamics (and of course typically 

S> A). In fact microscopic calculations of for the Pe-8 moleculeoLJ can be done (recall Fig. 5(b)), since the 
hyperfine couplings are essentially dipolar (the relevant nuclei include 120 protons, 8 Br nuclei, and 18 N nuclei). 
One finds that even the weak molecular dipolar fields strongly distort the hyperfine fields, mixing up the different 
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polarisation groups, and also-rijting a large value of k. The final value will obviously depend sensitively on any 
nuclear isotopic substitutionOLj. The subsequent r-«volution of the hole depends on sample shape; this has been 
studied theoretically using Monte Carlo simulationsE3. 

(ii) The short-time relaxation for the total magnetisation M{t) = J d^AI{£^,t) should have a " square- root" time 
form 

M(i7o,t) = Mo[l - (t/TQiHo))'^^] , (5.3) 

where Mq is the initial magnetization Iqij_with appropriate modifications for other protocols, such as a zero-field 
coofing followed by relaxation in a fielcHa), and Tq^{Ho) — c{£,o/ Ed)^"^ M{£, — —gfiBSHo,t = 0), where W^^ is 
the "density of states" of the distribution. The constant c is dimensionless, depending on sample shape, and can be 
evaluated analytically or numerically. Notice that the existence of the square root does not depend on sample shape, 
and indeed its raei i siate nce over fractional relaxations of ^ 0.1 is clearly demonstrated in Monte Carlo simulations for 



different shape£§E3^EI. 

This result implies that by varying Hq, one can measure M(^), by extracting Tq^{Hq) at successive values of Hq. 
If one knows M(^) (as one does in an annealed sample - it should be Gaussian) then one may then extract A from 
measiir eiTierits of Tq^ . The \/t dependence of M{t) has since been reported in quite a few experiments, on both Fe(8) 
crystalsOH'Ea and Mn(12) crystalsO (although the situation in Mn-12 is seriously complicated by " impurities"c3) . 
The Fe-8 experiments have produced reajiM-kable "maps" of M(^), and its time variationE3c3. The "hole-digging" 
has been found in both Fe-8 and Mn-12c3cJ, with an "intrinsic" short-time intrinsic linewidth which is roughly that 
expected from the hyperfine interactions, provided one takes account of the effect of internal fieldsE3, which make Ti 
very short (so that s-is effectively ~ \fN)- 

Wernsdorfer et aE3 alsp used this technique to extract the value of A for Fe-8. In a remarkable experiment, 
Wernsdorfer and SessoliE3 extended these measurements of A to include ji transverse applied field H^; as noted in 
section 2.F, the topological giant spin phase should vary with field -ffj_LZnIj, producing Aharonov-Bohm oscillations 
in A(_ffj^). These oscillations were found, both as oscillations in the relaxation rate Tq^{H±), and using a quite different 
AC absorption ("Landau-Zener") technique; these independent techniques agreed rather well in the measurement of 
A. Notice, incidentally, that neither method can properly measure A near its nodes (ie., where A — > 0); this is because 
of the distribution of internal transverse fields. Incidentally, we should strongly emphasize that these experiments 
(even the Aharonov-Bohm ones) do not demonstrate coherent tunneling- indeed they show exactly the opposite! This 
is because the experiments are inherently relaxational (this is why all rates are ~ |Ap, and not ^ A). Readers puzzled 
about how an Aharonov-Bohm effect can occur in a relaxation rate, are encouraged to think about other examples in 
physics where phase interference shows up in irreversible quantum phenomena. In the present case we may loosely 
define a "decoherence time" for the molecular spins (one should not push this too far, given the non-exponential 
nature of the central spin relaxation!), and one finds that t^Aq <C 1 (no coherent oscillations possible) but r^fio ^ 1 
(ie., coherence is maintained diping a single very rapid tunneling transition). 

In a very recent experimentEil, deliberate modification of the nuclear isotopes in an Fe-8 crystal has shown a 
dramatic modification of the bulk magnetic relaxation- this is the most direct evidence so far for the role of the 
nuclear spins in mediating the quant una, jxelaxat ion. It will be interesting to see a quantitative comparison with the 
calculated results for different isotopesE^CJ. 

(iii) Quantum Spin Glasses: The work on nanomagnets probes our understanding of collective phenomena in many 
magnetic systems, ranging from spin chains to quantum spin glasses. This latter exa mple is very closely related to 



the molecular nanomagnets just described, since the model Hamiltonian is just (5T), with the dipolar interaction 
replaced by a set {Jij} of " frustrating" interactions, usually long-ranged. A typical experimental example is provided 
by the disordered dipolar spin system LiIl0:rYi_2:F4 (where the Ho moments interact via dipolar interactions^, and 
the A are induced by a transverse field). Until now most theory has ignored the environment in this problem, but 
from our discussion above, this is clearbL,a mistake if one wishes to discuss dynamics. This is also the view expressed 
in recent papers of Cugliandolo et al.E3, who have included coupling to an oscillator bath environment. This is 
pres umably correct in metallic glasses, where the coupling to electrons dominates (compare the discussion following 



(2_^)), but in insulators the coupling to nuclear spins dominates (in Lii_2,Hoa;F4 this is clear- the hyperfine coupling 
even modifies the phase diagram!). An observation of strong hole-digging in M(^, t) in a spin glass would be consistent 
with this, since oscillator baths typically give most rapid relaxation away from resonance {of. discussion after eqtn. 



(3.45). Note also that one of the main theoretical questions concerning quantum spin glasses, viz., what happens 
when A ^ Jij, could be examined experimentally using molecular magnets like Fe-8, or the LiHoa,Yi_:cF4 system, 
in a strong transverse field, where one might reasonably expect to see coherent propagation of tunneling events from 
one site to another if A is sufficiently large. 
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B. Conductors and Superconductors 



In mesoscopic conductors the standard weak locahsation theoryO evaluates a " decoherence time" (T) in terms of 
the electron-electron scattering rate, which is itself strongly influenced at low T by elastic impurity scattering. This 
is essentially an oscillator bath problem- the relevant oscillators being diffusons, Cooperons, a*id phonons. Curiously, 
given the known importance of scattering off dynamic " 2-level" fluctuators in these systemsES, there has Jaeen little 
theory on the effect of these on r^, apart from the pioneering work of Altshuler and Spivak, and Feng et alo. This is 
of course a spin bath problem, with the spins representing defects, paramagnetic impurities, etc., in the environment. 
It will be interesting to see if the "saturation" in t^{T) reported at low TE3 may at least be partially explained by such 
scattering. The physics of this saturation depends on the more fundamental question of how decoherenre. behaves in 
the low-T limit in conductors, and has caused considerable debate in the recent theoretical literatureEJO. 

In the case of superconductors the situation is similar, in that almost all theoretical work on dissipation and other 
environmental effects has looked at the effects of electronic quasiparticle m(jidfis.p£|^otons, ie., delocalised modes 
which can be mapped directU Xo oscillator mod£a._The pioneering papersll2lll3cZI'Ell led to a massive subsequent 
literature, both experimentallirllj and,theoreticaE3'0 In most work on tunneling there is no question that theory 
and experiment correspond very wello. However the situation is more delicate for coherence, discussed below--, Iii 
spite of the extensive theory of spin impurity effects in superconductors we are aware of only 2 theoretical paper3l3't3 
(and no experiments) examining their effect on the flux dynamics (in particular tunneling) of SQUIDs. 

C. Coherence and decoherence; and "qubits" 

An understanding of decoherence mechanisms is central to the exploitation of mesoscopic systernsLia. quantum 
devices, as well as to general jqiiestions about how quantum mechanics applies on the large scale,Ea'L3, and the 
quantum measurement problemEa. It has become particularly important now that efforts are being made to construct 
"qubit" devices, with a view to making quantum computers. 

The last 15 years have seen a total transformation in how such questions are discussed- instead of vague analyses 
in terms of "measurements" by the environment, we now have precise and generally applicable models, which can 
be tested in many experiments. The basic issues are (i) whether phase coherence can be preserved in the reduced 
density matrix of the system of interest and (ii) what are the decoherence mechanisms destroying it. Here we briefly 
review studies of superconducting and magnetic systems, and then examine things from a more general theoretical 
standpoint. 

(i) DecoherexLCE. in superconductors: This has been discussed intensively ever since the theoretical predictions of 
Leggett et al.E^pa concerning "macroscopic quantum coherence" in SQUIDs, and subsequent proposals for experi- 
mental searchesEll. To date no experimental success has been reported (althoueh-tliere is good evidence for resonant 
"one passage" tunneling transitions between near degenerate levels in 2 wellaljCj). Almost all microscopic analy- 
ses of thi s i piyph liSgi have assumed environments of electronic excitations which can be mapped onto oscillator baths 
(see, eg.jEaEDea'ca). In our opinion, as discussed in section 4.B, the basic problem is simply that the main source 
of decoherence in most systems (including SQUIDs) at low T will not be any oscillator bath, but the spin bath of 
paramagnetic and nuclear spins. As discussed in section 2.G and 4.B, the low-energy scale of this spin bath means it 
will not usually have a big effect on SQUID tunneling, but its effect on macroscopic coherence or on superconducting 
qubits will be rather large. i— . 

Although a superconducting qubit has not yet been built, experiments may be getting rather closecJ. To see how 
big spin bath effects on coherence might be, let us recall that the effect of paramagnetic impurities is to create a 
Gaussian multiplet of spin bath states of width Eq ~ fJ^sBo^/ Npm in energy for a SQUID containing a total of Npm 
paramagnetic impurities interacting with the supercurrent, where Bo is the change in field on each paramagnetic 
impurity caused by the change in flux state of the SQUID. To see coherence it is necessary, from the discussion of 
section 4.B, that Aq » i?o, because the decoherence time coming from the^pin bath is ^ Ao/E^; this essentially 
sets a lower bound for Aq. As discussed in some detail in a recent paperEa, jthis turns out to be a rather stringent 
requirement on real SQUIDs; in fact in the experiments of the Lukens groupCJ one infers a value Eq ~ 0.4 K from 
their resonant linewidths. Obviously this value could be reduced a great deal by careful attention to the nuclear and 
paramgnetic spin impurity composition in the system (as well as to the sample geometry). 

(ii) Decoherence in Mamets: The most dramatic claims for the observation of macroscopic coherence have been 
made by Awschalom et al.E^, working on randomly oriented dilute ensembles of ferritin macromolecules (which order 
antiferromagnetically, but carry an excess moment of somewhat random size; the antiferromagnetic "Neel" moment is 
~ 23, 000 /is). In an effort to make the molecule size as uniform as possible, these authors filtered them magnetically. 
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They also artificially engineered molecules of smaller size. The essential result was the observation of an absorption 
peak at MHz frequencies, whose frequency varies approximately exponentially with the size of the molecules. This 
was interpreted as a signature of coherent tunneling between "up" and "down" states of the Neel vector. There have 
been widespread objections to this interpretation, both on theoretical and experimental groundalHj, and so far no 
other group has succeeded in confirming the experiments. Note that the Awschalom group saw similar resonances 
(also with an exponential dependence of resonant frequency on size) in large FeCos particles, but did not attribute 
this to tunneling 
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As we discussed in the previous sub-section, there is now very extensive evideruje, that, nanomagnetic molecules in 
macroscopically ordered crystals tunnel incoherently in the low-T quantum regim^3C3~E3. There is thus an apparent 
contradiction between the ferritin work and that done in the Mn-12 and Fe-8 systems (particularly since the ferritin 
molecules are much larger and certainly contain a lot of spin disorder). Thus in the very well-characterised 
molecular crystals used by the Florence and Grenoble groups, the parameter n characterising decoherence from the 
orthogonality blocking mechanism variesEaL^, even in an ideal sample, between k ~ 6 — 15 in zero applied field 
(depending on the annealing-dependent spreadj-H^D in the intermolecular dipolar fields), to k ~ 80 when ^ 0.2 T 
(where the first zero in A is supposed to occuicj). Recalling from sections 3.C and 4. A that coherence is practically 
eliminated unless k <C 1, it is hard to see how experiments on these particular molecules in low fields will stand much 
chance of seeing it. 

On the other hand it is clear that future experiments on single nanomagnets in the quantum regime might have 
a chance of seeing coherence iff one could raise to values ^ the hyperfine couplings cj^ (presumably using a 
large external transverse field), thereby making k <^ 1 and so removing decoherence (and also reducing the problem 
to a straightforward spin-boson model- see section 4.B). Another possibility, which could be realised in, eg., the 
IaROx^ i-x^ A system at high transverse fields, would be to see coherent propagation of spin flips (ie., spin waves) in 
a lattice of spins, by making A > Wd (here one could also make uok ^ kT , thereby freezing the nuclear dynamics!). 
Conceivably the same could be done in an Fe-8 crystal (now with the inequality A > Wd 3> LOk operating). We see 
no reason why such experiments could not be done in the next few years. 

(iii) Decoherence as T Q: Let us now consider the general question of decoherence effects at low T . Decoherence 
is often (particularly in conductors) characterised by a "decoherence time" r^, for the phase dynamics of the degree 
of freedom of interest. If and when is meaningful, it may b£,much shorter than the energy relaxation time te {cf. 
the example of a single oscillator coupled to an oscillator bathEj, or the examples of topological decoherence given for 
the spin bath in sections 3. A and 4. A above)). Coherence exists if t^Aq <C 1, where Aq is the characteristic frequency 
of the system's phase dynamics. 

Notice that what allows us to discuss this problem with any generality at all is the assumption, discussed in sections 
1 and 2, that a few canonical models describe the low-T behaviour of most nbysical systems. Extensive study of the 
relevapt jeanonical oscillator bath models (in particular, the spin-boson modelt3^lijl£3 and the "oscillator on oscillators" 
mode]ll3'Eil) show that with a power-law form J{uj) ~ w", decoherence disappears as T = for n > 1; for the Ohmic 
form J(w) = TTauj decoherence is finite at T = 0, but can be made small if a <g; 1. If the electro»ic-spectrum 
is gapped the Ohmic dissipation falls off exponentially in the low, 21 limit (thus for superconductorsc3c3 one has 
J(tj,r) - we-^scs/fcT fQj. ^ ^ 2Abcs, and for magnetic solitonsE^QH one has J{uj,T) - w(/cT/A„)e-^'"/'='^ for 
Lu < 3A„i, where Abcs and Am are the BCS and magnon gaps respectively). Thus, if one believes the oscillator 
bath models, coherence ought to be easily observable at low T; the condition r0(r)Ao ^ 1 is clearly satisfied for 
temperatures well below the gap energy. _ _ 

If we examine the canonical spin bath models we find a very different storytTu. Consider first the central spin 
model as T — > 0; we will go to such a low temperature that all of the spins in the spin bath order in the field of 
the central spin (ie., T IfiK in some cases), ie., all intrinsic fiuctuational dynamics of the bath is frozen out. 
Does decoherence disappear? No, because the mechanisms of topological decoherence (induced bath spin flip) and 
orthogonality blocking (precession of the bath spins in between central spin flips) still exist- the bath can still acquire 
dynamics from the central spin. We emphasize here that this physics cannot be described by an oscillator bath model. 
From the results in section 4 we see there is a residual constant decoherence as T — > 0, coming from the spin bath, 
which in an experiment would be signalled by a saturation of once oscillator bath effects had disappeared. The 
extent of this decoherence is characterised by the behaviour of PM=oit), (cf. section 4. A), and we saw that unless 
both K and A were ^ 1, decoherence was strong. In the case where all bath spins are polarised by the central spin, 
transitions are blocked anyway- there are no spins in the M = polarisation group! 

We conclude that for any system described by the central spin model (ie., where the central system reduces to 
a 2-level system at low energies), a general consequence of the coupling to a spin bath will be a loss of coherence, 
via either the topological decoherence or orthogonality blocking mechanisms, even in the T — > limit. The residual 
coherence (if any) will depend on the strength of the couplings to the spin bath, in a way discussed quantitatively in 
sections 4 and 5. 
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We may generalise these considerations to models in which a " particle" moves through a spin bath (the same model 
also describes a network of spins, or of mesoscopic superconductors, etc., coupled to a spin bath), and get the same 
result. Conaider, eg., a particle hopping from site to site on a I?-dimensional hypercubic lattice, whilst coupled to a 
spin bathE23, and described by a Hamiltonian 
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where < ij > sums over nearest neighbour sites. The couplings "w^:' and Vj^-' to the bath spin ak usually depend upon 
which "lattice rung" (i, j) the particle happens to be, because the coupling normally has a finite range and the spin 
(Jk has some position with respect to the lattice (one can also add a longitudinal coupling to the {i7|} on each site, 

as in eqtn. ( 2.30|) ). In the continuum limit ( |5.4| ) reduces to ( 2.11 ) and ( 2.12 ), after dropping the external field hk 
and the dependence of F^{P,Q) on Q. To isolate the decoherence effects let us assume Vkk' = 0, ie., we now study 
the analogue of the strong-coupling regime in section 4. A. Without loss of generality we may then concentrate on the 
"phase decoherence" regime, ie., on an effective Hamiltonian i?e// — X^m ^mH'jJ^^ , where 
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This is just a generalisation of (4.5) to the lattice; the dependence of the coupling on the lattice position is dropped 
because it is inessential to what follows. 

Coherence, if it exists, will appear in the function Pno{t) (the probability to start at site 0, and be at site n a time 
t later). Perfect coherence (ie., no bath) yields 
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for which {t) ~ l/(Aot)^ at long times; moreover, the 2nd moment {{n{t)'^) = X^n "^^so'C*) D{Aot)'^ at long 
times. Both are characteristic of coherent "ballistic" band motion. 

If we now go to the interacting case-one finds, by a straightforward generalisation of the calculations in sections 3 
and 4, some rather interesting resultsll23. Consider first a bath with all spin states equally populated. Then at long 
times ((n(t)^) ^ D{Aoty but PQo\t) (1/Aot), independent of D\ Mor cover, if we start with an initial Gaussian 
wave-packet of width An{t = 0) Ro, one finds Pgo^t) (l/RoAot) as t oo. These results show that the naive 
inference of ballistic propagation from the second moment result is wrong- in reality one has strongly anomalous 
diffusion (with an energy-dependent diffusion coefficient, demonstrated by the i?o-dependence of the results). In fact 



(0)/ 



the probability function Pfi^it) decays like Ro/{Aot)\i 
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at long times (ie., Aot ^ Ro) and distances \n\ ^ Ro, 



out to a "ballistic" distance l{t) ~ Aot, for dimension D > 1. Thus there is an advance "ballistic front" which decays 
in amplitude with increasing distance/time from the origin, but only as a power law; and it is followed by a much 
larger anomalously diffusive (and of course incoherent) contribution. For sufficiently long times (in fact for AAq^ ^ 1) 
this result is independent of A. Similar results apply for the orthogonality- blocked case as a function of k. 

Now let's take the limit T ^ 0, meaning that we allow the bath spins to order in the field of the particle. It is 
clear that if aj| 3> Aq for all spins, and provided there are no non-diagonal momentum couplings to the bath, we 
cannot get decoherence by the same mechanism as above, since there is only one state in the relevant polarisation 
group- the particle will then move freely without disturbing the spins in any way. However in any other case phase 
will be exchanged with the bath in the same way as above, with or without dissipation (which will certainly arise in 

the weak-coupling limit wj! < Aq). 

We therefore conclude that a finite decoherence in the T ^ limit is a generic consequence of the existence of spin 
bath environments. 

(iv) Qubits and Quantum Computation: What is the impact of these results on hopes for quantum computation? 
A Quantum computer is an information processing device which can be imagined as an assenjibly of 2-state qubits, 
these being none other than the spin-boson and Central spin models discussed in this articlai23'E£j. Such a computer 
has yet to be built, and papers and books on this topic tend to divide into 2 classes. The first simply ignores 
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decoherence (apart from occasionally referring to it as the main stumbling block preventing the construction of a 
quantum compute*!^, whereas the second regards decoherence as the crucial problem, and eithei^ries to treat it 



theoretically (eg.,liH3) or maintains that a quantum computer will never be built because of it (eg.,t2II) 



The basic problem here is the lack of serious theory on the effects of decoherence, starting from realistic models 
which can be tested quantitatively by experiment. The analysis of the Central Spin model in sections 3 and 4 can 
be used for a single qubit- it will be clear that the main task is to reduce diagonal couplings to the spin bath (and 
also any oscillator bath) as far as possible. However this is only the beginning of the problem- the operation of a 
quantum computer involves multi-qubit wave-function entanglement, and thus one wants to understand the behaviour 
of a decoherence time t^^^^ (^i, ^2, - -Cm), governing loss of M-spin phase correlations, in the presence of coupling to 
spin and oscillator baths. We are aware of no studies of this problem (or even recognition that it is a problem) m 
the literature. Studies of mutual coherence and decoherence in the problem of 2 spins coupled to an oscillator bathE2l 
give some feeling for what might happen, as do the studies of lattice systems. However we still have no answer to, 
for example, the question of how t^^^^ behaves for large M . If it ends up having a generic decrease ~ e~°*^ then the 
whole quantum computing enterprise will be in very serious trouble! 

It is clear that in the near future research on quantum computers will have to proceed on both practical designs 
(with serious attention given to the decoherence characteristics of the relevant materials) , and also on general studies 
of decoherence for systems of coupled qubits, themselves coupled to spin and oscillator baths. This promises to be 
one of the great challenges in condensed matter physics during the next decade. 
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VII. NOTE ADDED IN PROOF (FEB 22, 2000) 

A number of papers touching upon the present subject (particularly on decoherence) have appeared since this article 
was written. . .. . 

On the highly controversial question of zero-temperature decoherence in mesoscopic conductorsOEJ (which accord- 
ing to some authors throws the entire conventional theory of metals into doubt), i-a-large number of suggestions have 



appeared. A list of them appears in a short experimental review by Mohantyll23,p4if|hicli_concludes that only the 



original suggestion, of electronic coupling to zero-point fluctuations of the EM fieldEjaliilil, can explain the exper- 
imental jSatjUEation of r^. This conclusion is hotly disputed by other authors, both on experimental and theoretical 



grounds The idea that the decohercBce might be coming frona-^wo-level systems" in the sample (ie., from a 

spin bath) has been explored by Imry et al.tiS and Zawadowski et al.tHl, but rejected by Mohanty, mainly because it 
would imply very large low-frequency noise levels in the sample. In our opinion this requires further work- some of 
the decoherence mechanisms discussed in the present article do not appear in these papers, and would not necessarily 
show up in the low-frequency noise. In any case this controversy shows clearly how the theory of transport, including 
weak localisation theory, depends crucially on a correct understanding of decoherence mechanisms. 

The most interesting recent experimental progress concerning j-the mechanisms of decoherence may be work in 
magnetic systems. The most recent work by Wernsdorfer et al.Ed, which looked at the dependence of tunneling 
rates and "hole widths" in the resonant quantum relaxation of crystals of _Fe-8 molecules, has been supplemented by 



further work on the same systemEiS. Taken together these experiments give rather strong evidence that the tunneling 
is mediated and controlled by the nuclear spins in the system. No direct test has yet appeared of the theory of the 
decoherence coming from these same nuclear spins, however (which depends on calculati©ns-|Of the parameter k; cf. 



section 2.6, including Fig. 5, and sections 3.3, 4.1, 5.1, an(i-5,3, including Fig. 8, and refs.EfHIl]). Experiments looking 



for coherence in _Fe-8 have also been done at higher fieldsllij; however we are unable to see why the reported results 



give evidence either for or against coherence, since they simply show very broad and rather weak peaks in the ESR 
spectrum as a function of field. We emphasize a point here which has often been made, viz., that any demonstration of 
coherence requires direct observation of combinations of multi-time correlation functions- even a very sharp peak in, 
say, the AC absorption is not enough to demonstrate coherence. In this connection the reader is referred to Figs. 7-9 
in the present paper, which show peaks in x"('^) even when there is no coherence at all. Several reviews of the many 
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different experiments in tunneling nanomagnets have also recently appearedlli3~Lia. One topic not covered in these 
is the LiH0a;Yi_2;F4 system (section 5.1.3); some interesting new results-pn this compare the thermal and quantum 
annealing, showing the efficiency of the latter in quantum optimisationllia. p 
More general discussions of coherence and decoherence appearing recently include several theoretical reviewsE 
(although these do not really discuss spin bath environments) . There has slso. been interesting experimental work on 
systepis, other than superconductors and magnets; see, eg., Wiseman et al.E£3 for coupled quantum dots, and Myatt 
et al.llH3 for decoherence in quantum optical systems. 

APPENDIX A: DERIVATIONS FOR LIMITING CASES 

In this appendix we give the derivations of some key formulae in sections 3 and 4. Instanton methods are used to 
handle the spin environment. In Appendix A.l we give details of the fairly trivial calculations required to deal with 
averaging over bias, topologi cal p hase, and fluctuations in bias. Then in Appendix A. 2 we discuss the more lengthy 



derivation of the expression (3.21) involved in orthogonality blocking. 



1. Topological phase, bias, and bias fluctuation effects 

We wish to evaluate Pii(t), the probability for the central spin to return to an initial state | t) after tim e j . , in 
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T 

the presence of a static bias ^ and a noisy bias e(t). We use ssh physics. ubc.ca standard instanton techniquest2 
because they easily generalise to include the spin bath. 

We begin by ignoring the topological phase of the central spin since its effects are trivial to add. Then the amplitude 
for a 2-level system to flip in a time dt is iAodt, and the return probability is given by summing over even numbers 
of flips: 

A°i(i)-^|l + f:^^^^| = i[l + cos(2A„i)], (Al) 

Now consider the modification introduced by a longitudinal static bias f . In instanton language we begin with the 
return amplitude Aii{t,£,), and Laplace transform it; using the action e**^''* in bias ^, over time dt, we get 
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where Eo — [S,^ + ^l]^^^ ; then the standard answer for Pii(t, ^) is just: 



Pii(t,a=/. dp,dP2 ._^^^._^^. j:Y: 
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[pi - iP2 - ^0 ^fr'o \pt+e J \pi+e 



= [1 - ^(1 - C0S2M)] = [1 - |f sin^M] (A4) 
This representation avoids the difficulty in the usual representation coming from the square root E = + A2]i/2 

m 

the cosine. 

We now add a fluctuating bias e(i) to ^, with the correlation {[e{t) — e(t')]^) = A'^jt — <'| for short times. Noise 
averages are then given by the Gaussian average 

{F[eit)]) ^ /"p^^^^ ^j^jg-i/dti/dt2e(ti)x(ti-t2).(t2) 



1 j ds€.(^s^ 

with a noise correlator 2iir^^(ii —t2) = A'^dtij + |i2| — |ti — ^2!). Averaging over a phase function F(ti, i2) — e 
then gives the standard result 
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(A6) 



In applying this to (A2) we assume fast diffusion (see text), and thus only expand to ^ O(A^). This gives Pii{t) = 
e-t/^(C) with 
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where G(e) = (2/7rr^)^/^e~^^'^/^^^ is the probability e{t) takes value e (ie., it is the lineshape of the polarisation group). 

Adding the the topological phase $0 = ttS to these calculations changes the flip amplitude to iAo exp{±i$o}(it; one 
then sums over all paths with an even number of flips and over all combinations of ± (clockwise and counterclockwise) 
flips; thus (Al) becomes 
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with similar obvious modifications to include bias. The generalization to include the phase from the bath spins 
(topological decoherence) is now obvious for both imaginary and full complex Uk (see text, section 3. A). 



2. Orthogonality blocking effects 



We consider the situation described in section |^.B., where the "initial" and "final" fields on the bath spin are 
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'f^"'^-' and 7^^\ related by an angle (3k, which is assumed small, and is defined by cos2/3fe = —7^^'' ■ l)^' /\lk 
choose axes in spin space such that the initial and final spin bath wave-functions are related by 
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where tJk = e 



(compare eqtn. (3.19) 



In general the initial spin bath state will belong to some polarisation group Mo (not necessarily Mo = 0), where the 
polarisation is defined along some direction defined by the central system (for example, in a nanomagnetic problem, 
one could define it as the direction of initial orientation of the nanomagnetic spin). As explained in the text, during 
a central system transition energy conservation requires the polarisation to change from Mo to Mo — 2M (and back, 
for further transitions); for "pure" orthogonality blocking (ie., when no other terms are involved in the Hamiltonian) , 
Mo = M. In what follows we calculate the correlation function PMo,M{t), the central spin correlator defined under the 
restriction that the spin bath transitions are between subspaces defined by {V) = {J2k=i ^k) — and {V) — Mq — 2M 
subspaces, which are supposed to be in resonance. The statistical weight of states with Mo > N^^^ is negligible, so 
we will assume that Mq, M < N. 

We enforce the restriction to a polarisation group M using the projection operator 
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We can now write down an expression for the amplitude (not the probability!) A^j ]^.j{t) for the central spin t to stay 
in state | 1) during a time t: 
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where T2„ is 
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^From (All) we can now write the full correlation function Plj ^^{t) as 
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(Tj^r2„) . (A13) 



We now use the assumption that the (3k are small; mo re precisely we assume that the orthogonality exponent 
K, defined previously by e^** = J^cos/3fe (cf. eqtn (3.23)), can be approximated by the perturbative expansion 
K « i ^ 0^. This assumption makes it much easier to calculate the average in (A12). We consider first the problem 

with only one environmental spin ak, and calculate the average {T2^T2n)k in this case; since rJm^^n ^ product of 
operators acting separately on each a^, the average over all spins is also the product of single spin results. 

We only need consider processes with 0, 1, or 2 flips of the environmental spin, i.e., we expand in powers of Pk, 
and stop at Then it is clear that, if the initial state of <Jk is |tfc) 
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where the first term arises from the sequence [11 "f . . . 11], the second from the sequence [11 t • • • Tiii • ■ • with a 
flip when j — /; and so on. In the same way we find 
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to order (3f.. The sequence (|fe| {T2mV'^2n life) '^i^l have a similar expression, but with reversed signs coming from 
the 6*^^"^^ factors. 

We now observe that the state with polarisation Mo consi sts of [N + Mo) /I spins up and {N — Mo)/'2, spins down. 
Consequently, for each a^, we add f or i averages like ( A15), and then take the product 



N, 



k=l 



k' = l 



Substituting ([A15| ) into this expression we get 

(rLT'2„) =e'''°^^'"^'"^r«^)exp{-if,l/(e.,e;,M<,)} , 

where the "effective action" K'^H (^i, ^j,Mo) has two contributions K'^-f^ = Ki + K2: 

I'-i p'-i 



(A16) 



(A17) 



M ( 

Ki = 2k{1 -^)[{n + m) + ^(-1)''"' cos[2 ^ 6;] + ^ (-1)^''-^' cos[2 C^] 

^ /'>/ i—l p'>p j—p 

2m 2n 2n p—1 

-EE(-i)'"^-^[2E^^-Ee;]l 

p=l 1=1 i=l ] = l ^ 



(A18) 



34 



I'-l p'-l 



i'>; i=i p'>p j=p 

2m In p— 1 2n n 

- E E(-i)'^' ^^p[2* E - 2^ E ' (A19) 

We recall now that Mo < N^l'^ <^ N, which allows us to neglect the contribution due to K2 an d dro p the correction 
^ Mo/N to the coefficient k in Ki. Notice also that the phase factor in front of exp{— if'^^-'} in ( AI6 ) cancels exactly 
the phase proportional to Mq in the formula (A13) for PMo,M{t). Thus, quite surprisingly, we find the correlation 
function to be independent of Mg in this limit: 

pM(t) = E E (Z!Lv n n / f / ^ -p {2^M(6„_, + 6„-3 + . . . + 6) - e;)} , (A20) 

n=Om=0 ^ ''^ '' i=lj = l"' 

We can render this expression more useful by changing variables; first we consider the whole sequence = 
(^1, . . . , ^2rt, —^1, ■ ■ ■ , — C2m) together, and then define new angular variables 



2(n+m) 



E 2ea'+7ra, (A21) 



so that now 



exp|zAfE(-l)"+'Xa-2K[(n + m)+ E cos(xa - Xc)] | 

^ a a'>a ^ 



(A22) 



Thus we have mapped our problem onto the partition function of a rather peculiar system of spins, interacting via 
infinite range forces, with interaction strength 2k. 

To deal with this partition function , we define "pseudo-spins" Sa — (cos Xa, sinxa) and 5, such that 

2(n+m) 

5= E Ecos(^"-^"') ' (^23) 

OL—l a',a 

We can think of Sa as rotating in our fictitious angular space defined by the projection operator ( [AlOj ). Now consider 
the term G(5) in (|A22[ ) defined by 

/2(«+m) \ 

Gi^^) = n y ^e^^'^-'r-\. exp { - « E '^"^(Xo -Xa')} 



a = l 
'2{n+m) 

2n 



/ dx^^^Mi-lr^\^ I ^ (A24) 



Q = l 



This is easily calculated, viz.. 



2(n+m) ^ 

aXa »M(-l)° + iy„ 



a=l 



GiS) ^ I dSe''^^^- n /^e^'^(-i)°"^^"5(5-E^^a) 

a 

fSir J,, \ 2(n+in) 

/O 



2tt I l./n 27r 



^ fdzze~^-'/*-jl^+^\z), (A25) 
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where Jm(A) is the Mth-order Bessel function. Using 

E^2(n+m),2s _ Sgfi + 2^^ 
(2TO)!(2n)! ~ 2(2s)! 



(A26) 



to reorganize the sum over n and m in (A22) and changing the integration variable z — > 2xy/K, we then find 

[2iiA„($)JAf(2xy^)]2A 



poo 1 / 

PMit) = 2 dxx e-^' - f + V ■ 
■^0 V .=0 



poo poo 

PmH)^ / da;a;e-"='(f + cos[2Ao(«')JM(2a;\/;^)t]) EE 2 / dxx e'^'V/"^ (i, Aa/ (x)) . 



(A27) 



A A/ (a;) = A„(4>)JAf (2x7;^) 



(A28) 



Here we come to the crucial point in our derivation. Eq.( A27 ) gives the final answer as a superposition of non- 
interacting correlati on fu nctions for effective tunneling rates AA/(a;) with the proper weighting. For M = this is the 
form quoted in Eq.( |3.21 ) of the text. 

It is worth noti ng that non-zero M enters this calculation as the overall phase factor which we can follow from 
( A13| ) up to (A25), where we finally integrate over {xa} to produce the Bessel function of order M. This observation 
allows one to generalise any calculation done for M = to finite M by simply replacing Jg — > Jm in the final answer 
- the prescription which we use in other Appendices. 



APPENDIX B: DERIVATIONS FOR THE GENERIC CASE 



We outline here the derivations for section 4, in which topological decoherence, degeneracy blocking, and orthog- 
onality blocking are all simultaneously incorporated (the average over bias fluctuations being essentially trivial- see 
Appendix A.I). We have demonstrated in section 3 and Appendix A how each different term in the effective Hamil- 
tonian (2.16) influence the central spin dynamics. From these limiting cases we learned that static (or diagonal) 
terms in the Hamiltonian can be partly absorbed into a redefinition of the transition amplitude between states with 
equal initial and final energies. If we now deal with th e fu ll central spin Hamiltonian, we can still write the instanton 
expansion in central spin transitions in the form (see (A4)): 



PM{t) 



dpidp2 



= (Pl+P2)t 



oo oo 



{pi - ie){p2 



ie) ^ ^ 

' n=0 m=0 



(-»A,) 
pI -t- e 



-iA, 



Bnm{M) , 



(Bl) 



{9i=±} 



(B2) 



where the sum over {gi = ±} with 1 < Z < 2(n + m) describes all possible clockwise and anticlockwise transitions, 
and the operator product is defined as 



f2n{M,gi) = UM{gi)ult{g2)UM{g^) ■ . . uUg: 



M\g2n) , 



(B3) 



Unig) = fine 



M ■ 



(B4) 



Here IIm as before projects on the polarisation state AA^ = M, and /3k describes the mismatch between the initial 
and final nuclear spin states. If all the couplings were equal {ujk — ujo ^ Aq) then the above set of equations would 
be the complete solution of the M polarisation group dynamics. One may then further average over different grains 
in the ensemble by summing over different polarisation groups with proper weigthing. If there is a small spread in the 
nuclear hyperfine couplings, it will produce an internal bias field acting on the grain, as described in the text (section 
|^.C. The final answer for the ensemble of grains is obtained then by averaging ( pl| ) over the bias field. This bias 
field is due to all environemntal spins interacting with S, and hardly changes when a few nuclei (of order ~ A) flip 
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with S. For this reason there is no back influence of the induced nuclear spin flips on the bias field, at least during 
the time scale set by the damping of coherent oscillations (when many environemntal spins are flipped the coherence 
is obviously already lost). 

The crucial observ atio n is that if the sum over the clockwise and anticlockwise trajectories and the average of the 
operator product in (p3) can be presented as some weighted average and/or sum, of form 



(Af) = I dxidx2 . ■ -dxa ^ ZM{xi,...,Xa;ki,...,kb)R^J;P^™'\xi,...,Xa;ki,. 

^ 7, J, J. 



(B5) 



kik2 ■ ■ -i^b 



with flxed integer values a and 6, then the problem may be considered as solved because the instanton summation 
then reduces to that of a coherent (non-interacting) dynamics with the renormalized tunneling amplitude 



Ao — > AA/(a;i, . .. ,Xa; fci, . 
and the final answer acquires a form 



,kb) = AoRm{xi, . . . ,Xa;ki, . . .,kb) 



PM{t, e) = / dxidx2 ...dxa ^ Zm{ 



xi,...,Xa;ki,..., kb)P[i^^{t, e, AA/(a;i, ...,Xa\k 



,kb) 



(B6) 



(B7) 



ki k2...kh 



where P^^^^ {t,e, Am) is described by Eg. ( |A4| ). 

We have already seen that Eq. (B5) is indeed valid for the cases of pure topological decoherence and pure orthog- 
onality blocking- we now prove that it also holds when we combine the effects of topological decoherence with the 
projection on a given polarisation state. Here we evaluate the M — contribution; the result for PM{t) then follows 
from the generalisation explained at the end of the previous Appendix. 

Introducing as before the spectral representation for the projection operator [see Eq. ( |A10| )] we write 



{ai=±} 



2{n+m) 



(B8) 



exp {-KUilai}, {Cp})} = {e'^'^e"' ^A-^'' e'^2V ^^a2 g*«2(„+™)^g*92(„+™) ^A-ff,^ (gg) 

With the usual assumption t hat t he individual ak are small (but not necessarily A), the "effective action" K^(^ 
can be written as (compare Eq.( A22| )) 

2{n+m) 2{n+m) 

KUi{Cp})^X' W+A Y cos(xp-Xp').9pV . (BIO) 

p',p p',p 

which generalizes from orthogonality blocking; the xp ^re defined as in ( |A2ll ), and 

N N 



k=l k=l 



as before. We use the same trick of introducing "pseudo- vectors" Sp — (cosxp, sinxp), and = ^X]p=i^™^ .9p'^p) = 
J2p' p9p9p' <^os(xp — Xp'); to decouple integrals over the new variables Xp = Xp + '^9pl'^- After some lengthy, but 
straightforward algebra we get P\\{t) in the form 

Po(i) = 2 / te-^' E Fx'{mY^'^^'^-'^^PfXt,Ao{^,x)) (B12) 

where Aoi^p^x) — 2Ao cos (pJo(2a;VA) as before. 

The case of complex (even assuming uj-^ = in the effective Hamiltonian) is more subtle technically, but goi 
through in exactly the same way. Here we just outline the key steps; a more detailed derivation may be found ir 
The effective action now has the form 
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2(n+m) 



p.p'=i ^ 

+ cos(xp - Xp')[(A - A') - {-IY+P\^i - ry') - *(7 - l'){{-lY + (-1)"')] \ , (B13) 



where the constants are defined by: 

A 



AT , N 



2 ^ ^ ' 2 

k=l fc=l 



JV JV 

fc=i fc=i 



1 ^ 1 ^ 

7 = ;t E • ; = 9 E ^k^knlvl ; (B16) 



2 ^ ^ ' ' 2 

fc=l /c=l 



As before we change variables according to Xp = Xp + when gp — 1, to introduce odd and even spin fields 

2(ri+m)-l 2{n+m) 
p=odd p=even 



which are used to decouple integrations over Xp with the final goal to get the answer in the form of Eq. (|B5|). This 
indeed can be done, and the final answer reads 



27r 7 27r 

mi— — CX2 m2— — OO 



X Z{(pi, (p2, xi, X2,mi,m2)Pii' [t, Ao{xi, X2, (fii, (P2)] , (B18) 



and (B18) has an obvious generalisation to include the bias integration J de. The weight is given by 



2^ _ g2i[mi(<i> — — m2ip2+4mim27']g4(77'm2— A'rrij) 

X1X2 J f{a + b)xiX2\ { {a — b + 2ic)x\ + [a — b — 2ic)x'2 



^ Wb^)'' \ Kab~c^) ) 1 16{ab-c-) / ' ^^^'^ 

and the renormalized tunneling splitting equals 

Al{xi,X2,fi,f2) = 4A^cos((^i +(^2)cos(v3i - tp2)Jo{xi)Jo{x2) , (B20) 

where 

a = A-A'; b = t] - r]' ; = 7-7'. (B21) 
FIGURE CAPTIONS 

Figure 1 The flow of a class of effective Hamiltonians describing a central system coupled to a background 
environment, in coupling constant space, as the UV cutoff in the joint Hilbert space is reduced from Ec to Qq- Here 
we show flow to a fixed point FP, in a simplified 2-dimensional space of couplings ai, 012, but one may also have fixed 
lines or more complex topologies. 

Figure 2 A typical path for a 2-level central system (solid line) coupled to environmental modes (wavy lines) as 
a function of time, showing the couplings which exist in both the spin-boson and central spin models. We show both 
diagonal couplings D to and non-diagonal couplings ND to t± (in the central spin model these are strong enough 
to lead to multiple excitation of environmental modes). 
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Figure 3 Definition of tfie longitudinal and transverse parts of the diagonal coupling to a bath spin in the Central 
Spin Hamiltonian, in terms of the initial and final fields 7^^-' and 7^^-' acting on this spin- this also defines the angle 
(3, and the mutually perpendicular unit vectors I and rh (see text). 

Figure 4 Classifying the states of the Central spin Hamiltonian. Each level of f is associated with a 2^-fold 
multiplet of bath states (Fig 6(a)). These are classified into polarisation groups {M} (where M is the total polarisation 
along z), separated by energy cuo and with width Fm; Fig 6(b) shows the density of states Gm{0 of the separate 
groups, and Fig. 6(c) their sum W{S,). We show W{£,) for 2 different values of the parameter ^ — Tm/^o', in realistic 
cases » 1 (ie., the polarisation groups strongly overlap), and W{£,) is Gaussian. Longitudinal transitions between 
2 different polarisation groups Mi and M2 go at a rate Tj"^; transitions within a polarisation group at a rate Tj"^. 

Figure 5 Example of the application of the Central Spin model to a magnetic macromolecule (the Fe-8 molecule, 
further descibed in section 5). In (a) we show the effective tunneling matrix element |Ae//| = |Acos(7r5'-|-«/3ono.Ho|), 
for this easy axis/easy plane nanomagnet in the presence of a field Ho = xHx in the x-direction (transverse to 
the easy axis), assuming an angle ip between x and the magnetic "hard axis" (perpendicular to the easy plane). 
Aharonov-Bohm oscillations appear when Lp is small, so that the action of the 2 relevant paths on the spin sphere have 
similar magnitudes^ but almost opposite phase. For larger iy9, one path dominates over the other and oscillations are 
suppressed. In (b) we show a histogram of the w| for this system- the main figure shows the protons and the lower 
inset the N and O contributions. The upper inset in (b) shows the variation of Eg and with (the parameter 
is discussed in sections 4 and 5). These figures are adapted from Ref. [69]. 

Figure 6 Behaviour of Pii{t) in the case of pure topological decoherence. We show Pii{t) — 1/2 for intermediate 
coupling, for which -Pii(i) takes the "universal form" discussed in the text. 

Figure 7 The effect of relaxation on a statistical ensemble of central spins, each interacting with a spin bath. 
In (a) we assume that A = 0, k = 5 and N = 1000, and show the normalised time dependence of 3 different 
contributions Pnit) to the total relaxation function Pii{t); they sum to give a roughly logarithmic time dependence 
for the total function Pii{t). The small M contributions relax quickly (up to M ~ k), so the effect on an initial 
ensemble distributed over bias ^, at short times, is to dig a hole around zero bias, of width ~ kcJq. In (b) we show the 
spectral absorption function x"(^) foi' = 2, dividing this into the M — contribution and the contributions from 
higher M groups (which relax more slowly and thus peak at lower uj). 

Figure 8 The spectral absorption function xXf=o('^) ^'^^ several values of k, for an ensemble of central spins in the 
M — polarisation group, in the case where orthogonality blocking dominates, and degeneracy blocking effects (ie., 
a bias average) are also incorporated. Contributions from higher polarisation groups M ^ are not shown; they are 
spread over a range ~ Mcoo, up to ^ kuJo- Contributions from groups with M > k are negligible. 

Figure 9 Graphs of Xm=o('^) "projected topological decoherence" (ie., including a bias average over an ensemble 
in which topological decoherence dominates), for several different values of the parameter (A — A'). Contributions 
from higher polarisation groups, which are spread over an energy range ~ Awo, are not shown. 
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